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ABSTRACT 

A Chandra X-Ray Observatory ACIS-S imaging observation is used to study the population of X- 
ray sources in the nearby Sab galaxy M81 (NGC 3031). A total of 177 sources are detected with 
124 located within the D25 isophote to a limiting X-ray luminosity of ^3 x 10 36 ergs s~ x . Source 
positions, count rates, luminosities in the 0.3 - 8.0 keV band, limiting optical magnitudes, and potential 
counterpart identifications are tabulated. Spectral and timing analysis of the 36 brightest sources are 
reported including the low-luminosity active galactic nucleus, SN 1993J, and the Einstein-discovered 
ultra-luminous X-ray source X6. The nucleus accounts for ~86%, or 5xl0 40 ergs s _1 , of the total X-ray 
emission from M81. Its spectrum is well-fit by an absorbed power law with photon index 1.98±0.08 
consistent with previous observations (average index 1.9). SN 1993 J has softened and faded since its 
discovery. At an age of 2594 days, SN 1993J displayed a complex thermal spectrum from a reverse shock 
rich in Fe L and highly-ionized Mg, Si, and S but lacking O. A hard X-ray component, emitted by a 
forward shock, is also present. X6 is spatially-coincident with a stellar object with optical brightness 
and colors consistent with an 09 - Bl main sequence star. It is also coincident with a weak radio source 
with a flux density of ^95 /iJy at A = 3.6 cm. The continuum-dominated X-ray spectrum of X6 is most 
closely reproduced by a blackbody disk model suggesting the X-ray source is an ~18 M Q object accreting 
at nearly its Eddington limit. 

The non-nuclear point source population of M81 accounts for 88% of the non-nuclear X-ray luminosity 
of 8.1 x 10 39 ergs s _1 . The remaining (unresolved) X-ray emission is confined within ^2 kpc of the galactic 
center. The spatial distribution of this emission and of the resolved X-ray bulge sources closely follows 
that of the bulge optical light. In particular, there is no evidence for an X-ray signature accompanying 
the filamentary Ha or excess UV emission seen in the central ^$1.0 kpc of the galaxy. The shape 
of the luminosity function of the bulge sources is a power law with a break at ~ 4 x 10 37 ergs s -1 ; 
suggesting the presence of an aging (~400 Myr) population of low-mass X-ray binaries. Extrapolating 
this luminosity function to lower luminosities accounts for only ~10% of the unresolved X-ray emission. 
Spectroscopically, the unresolved emission can be represented as a combination of soft, fcT~0.3 keV, 
optically-thin plasma emission and of a T — 1.6 power law. The unresolved bulge X-ray emission is 
therefore most likely a combination of hot gas and of one or more large and distinct populations of low- 
luminosity X-ray sources confined in the gravitational potential and tracing the old population of bulge 
stars. The distribution of disk sources shows a remarkably strong correlation with spiral arms with the 
brightest disk sources located closest to spiral arms. The luminosity function of sources near the spiral 
arms is a pure power law (slope —0.48 ± 0.03) while that of sources further away exhibits a break or 
cut-off in the power law distribution with no high-luminosity members. This is interpreted as a natural 
consequence of the passage of spiral density waves that leave the brightest (when averaged over their 
lifetimes) and shortest-lived X-ray sources immediately downstream of the spiral arms. Consistent with 
model predictions, we conclude that the shapes of the X-ray luminosity functions of the different galactic 
components of M81 are most likely governed by the birth rates and lifespans of their constituent X-ray 
source populations and that the luminosity functions can be used as a measure of the star formation 
histories of their environments. 

Subject headings: galaxies: individual (M81) — X-rays: galaxies — X-rays: binaries — X-rays: stars — 
supernovae: individual (SN 1993J) 



1. INTRODUCTION 

Systematic investigations of the X-ray properties of nor- 
mal galaxies began in earnest with the Einstein observa- 
tory over two decades ago. The picture that emerged for 
spiral galaxies (see the early reviews by Long & van Spey- 
broeck 1983; Helfand 1984; and Fabbiano 1989) is that 



the bulk of the X-ray emission takes place in two distinct 
physical environments: the star-forming disks of late-type 
spiral and irregular galaxies and among the old stellar pop- 
ulation in dense globular clusters and compact bulges at 
the centers of early- type spiral galaxies. In addition to 
these trends along the Hubble sequence, variations were 
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sometimes found among spiral galaxies of similar mor- 
phological type suggesting a dependence on star forma- 
tion histories. In particular, the brightest X-ray emissions 
are associated with starbursts in merging and interacting 
galaxies (David, Jones, & Forman 1992). Thus, by tracing 
the endpoints of stellar evolution, the X-ray source popu- 
lations of external galaxies provide important clues to the 
physical nature and evolutionary history of their hosts. 

The contemporary view for spiral galaxies is rapidly be- 
ing refined following the launch of the Chandra and XMM- 
Newton X-ray Observatories. Moderately-deep images re- 
veal point sources to limiting X-ray luminosities of order 
10 37 ergs s _1 in galaxies out to Virgo cluster distances. 
While this samples only the high luminosity end of the 
distribution of X-ray sources, of order 100 sources are rou- 
tinely detected in normal galaxies similar to our own. Reli- 
able spectral analysis is usually limited to an even smaller 
subset of the brightest individual sources. Nevertheless, 
using probabilistic methods, the observed sample of X-ray 
sources can help us understand current-epoch galaxy evo- 
lution in its broader context. 

A formal expression of the relationship between the star 
formation history of local galaxies and their observed X- 
ray source populations has recently been put forth by Wu 
(2001; see also Wu et al. 2002a,b; Kilgard et al. 2002; Dal- 
ton & Sarazin 1995). There it was shown that the basic 
shape of the observed X-ray luminosity function is gov- 
erned simply by the birth and death rates of the source 
population under the assumption that the more luminous 
X-ray sources are shorter-lived. Thus, in the absence of on- 
going star formation, the luminosity function will develop 
a cutoff at high luminosity that evolves toward lower lu- 
minosity. Conversely, if the population of X-ray sources 
is replenished through star formation processes, such as is 
found in spiral arms, then a power law shaped luminosity 
function can be sustained. 

Certain complications arise when applying this basic in- 
terpretation to X-ray populations in individual galaxies 
(Wu et al. 2002a, b). Among these are the presence of dif- 
ferent classes of X-ray sources, such as supernova remnants 
and accreting compact objects, which evolve on differing 
timescales; alternative source-formation mechanisms un- 
corrected with stellar evolution such as binary captures 
in globular clusters; and non-steady or luminosity-limited 
emission characteristics such as those associated with X- 
ray transients and novae and in Eddington-limited neu- 
tron star binaries, respectively. Therefore, only when spe- 
cific counterparts to individual X-ray sources are identified 
and their multiwavclcngth properties assessed can the full 
power of the hypothesis of Wu et al. (2002a,b) be applied 
to address the nature and evolution of X-ray sources in 
different environments. 

The nearby Sab galaxy M81 (NGC 3031) is ideal for 
such a study in that it contains both a strong two-arm 
grand-design spiral pattern and a well-defined circumnu- 
clear bulge. The distance to M81, 3.6 Mpc, has been 
well-established from Cepheid measurements (Freeman et 
al. 1994) which are in good agreement with other dis- 
tance estimates (Ferrarese et al. 2000). Populations of 
several classes of objects in M81 have been investigated 
and catalogued including globular clusters (Perelmuter & 
Racine 1995; Chandar, Ford, & Tsvetanov 2001), super- 
nova remnants (Matonick & Fesen 1997), H II regions 



(Hodge & Kennicutt 1983; Petit, Sivan & Karachentsev 

1988) , and stars and star clusters (Zickgraf & Humphreys 
1991; Ivanova 1992; Sholukhova et al. 1998). In addition, 
the plane of the galaxy is oriented 32° from face-on al- 
lowing detailed mapping of the velocity field (Goad 1976; 
Rots & Shane 1975; Adler & Westphal 1996) for dynami- 
cal studies and testing spiral density wave models (Visser 
1980; Roberts & Hausman 1984). 

The center of M81 contains a compact radio core (Bar- 
tel et al. 1982) surrounded by a region of enhanced far- 
infrared (Rice 1993; Davidge & Courteau 1999), Ha (De- 
vereux, Jacoby & Ciardullo 1995), and ultraviolet (Hill et 
al. 1992; Reichen et al. 1994) emission extending to ~ 50" 
(~900 parsecs). This emission probably comes from an old 
population of hot, low-mass stars rather than from young 
massive stars (O'Connell et al. 1992; Devereux, Ford & 
Jacoby 1997). 

In contrast to the bulge, HI velocity contours show a 
sharp discontinuity beyond the bulge identified as a spiral 
velocity shock (Visser 1980). Downstream of this shock are 
regions of star formation in the spiral arms. The distribu- 
tions of these components are consistent (Kaufman et al. 

1989) with density wave models predicting a broad spiral 
density enhancement (e.g., Roberts & Hausman 1984). 

Beyond the visible disk of M81 is an envelope of neutral 
hydrogen (Roberts 1972) enclosing M81 and nearby group 
members. A bridge of gas, a relic of tidal interaction (Cot- 
trell 1977), connects M81 and the starburst galaxy M82. 

The hypothesis of Wu (2001) and Wu et al. (2002a,b) 
was motivated in large part by the initial results from our 
Chandra observation of M81 presented in Tennant et al. 
(2001). There it was shown that the X-ray luminosity 
function of the population of bulge sources displays a break 
at ~4xl0 37 ergs s _1 similar to that observed in M31 (e.g., 
Shirey et al. 2001). The X-ray luminosity function of the 
disk sources, on the other hand, follows a single power law 
slope over three decades in flux. This is what is expected 
if an impulsive episode of star formation occurred in the 
bulge in the past, ostensibly during an encounter between 
M81 and one of its companion galaxies, while continuous 
star formation in the disk is being driven by the passage 
of spiral density waves. 

Here we build upon the earlier work of Tennant et al. 
(2001). After presenting detailed information on the in- 
dividual X-ray sources in §|3| and in-depth analysis of the 
brightest objects in §[|, the properties of the bulge (§|J) and 
disk (§||) regions are addressed separately then discussed 
(§0) within the common framework of galaxy evolution. 

2. OBSERVATIONS 

The primary X-ray dataset is a 49926 second observa- 
tion of M81 obtained on 2000 May 7 with the Chandra Ad- 
vanced CCD Imaging Spectrometer (ACIS) spectroscopy 
array operating in imaging mode. Unless otherwise noted, 
references to X-ray data will refer to this dataset. The 
X-ray data was reprocessed by the Chandra X-ray Center 
(CXC) on 2001 January 4. This reprocessed data is used 
in this work. There are no significant differences between 
the reprocessed data and the originally-distributed data 
analyzed by Tennant et al. (2001). The observation was 
taken in faint timed exposure mode at 3.241 s-frame -1 at 
a focal plane temperature of — 120°C. Standard CXC pro- 
cessing has applied aspect corrections and compensated 
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for spacecraft dither. 

The primary target, SN 1993J, was located near the 
nominal aimpoint on the back-illuminated (BI) device S3. 
The nucleus of M81 lies 2/79 from SN 1993J towards the 
center of S3 in this observation. Accurate positions of 
these two objects and two GO stars located on device S2 
were used to identify any offset and to determine absolute 
locations of the remaining Chandra sources as well as ob- 
jects in other X-ray images and those obtained at other 
wavelengths. Table 1 shows that the positions are accu- 
rate to within 1.4". No offset correction was applied to 
the Chandra X-ray positions. 

TABLE 1 

M81 ASTROMETRY 



Object 




Chandra 


Position 


Catalogued Position 






RA 


DEC 


RA 




DEC 


SN 1993J 


9 


55 24.77 


69 1 13.4 


9 55 24.77 


69 


1 13.7 1 


Nucleus 


9 


55 33.19 


69 3 55.1 


9 55 33.17 


69 


3 55.1 2 


PPM 17242 


9 


55 1.00 


68 56 22.1 


9 55 1.00 


68 


56 22. 2 3 


PPM 17243 


9 


55 2.57 


68 56 21.2 


9 55 2.76 


68 


56 22. 1 3 



References.-(I) Marcaide et al. 1993; (2) Ma et al. 1998 
(3) Positions and Proper Motions catalog 

A charge transfer inefficiency (CTI) corrector algorithm 
(Townsley et al. 2000) was then applied to the Level 1 
event list to partially correct for the charge loss and charge 
smearing effects of CTI in the ACIS detectors. Matching 
response matrices were also provided by L. Townsley. Af- 
ter correction, a single response matrix is adequate for 
sources on the S3 device since the spectral resolution does 
not exhibit a strong spatial dependence. This is particu- 
larly advantageous for analysis of the unresolved emission 
extending over the ~8' diameter bulge of M81 which is 
located entirely within the S3 device. 

The corrector algorithm was applied to all but the front- 
illuminated (FI) device S4 because no correction was avail- 
able for that device. Instead, the Level 2 event list was 
used and the destreak algorithm 5 was applied to remove 
charge randomly deposited along pixel rows during read 
out. 

The entire dataset was then cleaned of bad pixels and 
columns and the standard grade set and events in pulse 
invariant (PI) channels corresponding to ~ 0.2 to 8.0 keV 
were selected for source detection. The range 0.3 to 8.0 
keV is used for spectral analysis. 

No periods of high particle background occurred dur- 
ing the observation. The 0.3 - 8.0 keV background in the 
BI device S3 is ^0.04 cts pixel -1 and that in the FI de- 
vices is ~0.01 cts pixel -1 for the observation. Separate 
background spectra were extracted from large source-free 
regions of each device. The background spectrum for S3 
was chosen far from the nucleus because excess X-ray emis- 
sion was detected near the nucleus (§S). The background 
spectra appear similar to deep quiescent blank sky compi- 
lations available from the CXC. 

The detector viewing area covers 57% of the optical ex- 
tent of the galaxy, defined as the ellipse of major diame- 
ter 26. '9 corresponding to the D25 isophote as tabulated 
in de Vaucouleurs et al. (1991), oriented at position an- 
gle 149°, and with major-to-minor axis ratio 1.94:1 corre- 
sponding to the 58° inclination angle of M81. This area 



includes all of the S3 device, approximately half of each 
of the S2 and S4 devices, and the outer corner of 13. The 
data from each device is analyzed independently owing 
to differing energy resolutions, low-energy responses, and 
background signals. 

In addition to this primary dataset, a 2.4-ks ACIS-S im- 
age taken 2000 Mar 21 and numerous ROSAT PSPC and 
HRI datasets were used to construct long-term light curves 
of the brightest sources (see Immler & Wang 2001 for their 
analysis of the ROSAT observations). 

3. THE DISCRETE X-RAY SOURCE POPULATION 

Table 2 lists the 177 discrete X-ray sources detected 
in the primary Chandra observation. The table lists the 
source positions (in order of increasing right ascension), 
the aperture-corrected number of source counts, the signal- 
to-noise ratio for the count rate, apparent visual magni- 
tudes derived from either Hubble WFPC2 images or from 
the Perelmutter & Racine (1995) catalogue of bright ob- 
jects, ACIS CCD device identification, global environment 
(where b denotes bulge, d denotes disk, and D25 denotes 
source is outside the D25 isophote), and the unabsorbed 
luminosity in the 0.3 - 8.0 keV energy range. The ta- 
ble also lists corresponding X-ray detections from Einstein 
(denoted "X", Fabbiano 1988) and ROSAT ("P" and "H" 
denoting PSPC- and HRI-identified sources, respectively, 
Immler & Wang 2001), and potential source type based 
on spatial correlations with catalogued objects or on other 
information. Explanations of the quantitative entries are 
given in the following subsections. References to Chan- 
dra X-ray source numbers in this work refer to the source 
numbering adopted for Table 2. 

3.1. X-ray Source Detection 

A source- finding method was used that assumes a source 
is located at a given position and compares the distribu- 
tion of detected events to a known point spread function 
(PSF). The algorithm first calculates the fraction of the 
PSF within each pixel within a detection region. Then, 
using the PSF fraction as the independent variable, it cal- 
culates an unweighted least squares fit of a straight line to 
the counts detected in the pixels in the region. If a source 
is present, then the slope of the line will be positive and 
will represent the total number of counts from the source 
(integrated over the PSF). The line intercept is the back- 
ground per pixel. A key value is the uncertainty in the 
slope and hence the number of source counts. The uncer- 
tainty is determined by applying the standard propagation 
of errors directly to the sums. 

The algorithm then calculates the estimated source 
counts (slope) and error at every pixel in the image. The 
estimated source counts divided by the uncertainty is the 
signal to noise ratio (S/N). If the S/N exceeds some thresh- 
old then there is a source in the neighborhood. However, 
since both source counts and error increase near a source, 
the S/N is not as sharply peaked as either component in- 
dividually. To best separate sources in confused regions, 
a source is defined to be a peak in the estimated source 
count spatial distribution and must also exceed a minimum 
S/N threshold. The threshold S/N is best defined by con- 
structing the S/N histogram. This histogram is roughly a 



5 bttp: / /asc. harvard. edu/ciao2.1/downloads/scripts/dcstrcak.p: 
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Gaussian core centered near zero (due to noise) with ex- 
ponential wings (due to sources). By fitting the core to a 
Gaussian, the S/N value for which the Gaussian will con- 
tribute less than one source in the field can be estimated. 
This threshold is as low as 2.5 for a source on axis and 
as high as 3.0 for a source far off-axis. For this paper a 
constant value of 2.8 is used as the threshold. 

The PSF used in the search can have any shape and 
the source-finding algorithm allows either a mathemati- 
cal model PSF or a high-fidelity simulated PSF available 
from the CXC PSF library 6 . A circular Gaussian approx- 
imation to the PSF does a good job of locating sources, 
in comparison to simulated PSFs, assuming the width of 
the Gaussian increased quadratically off-axis so that the 
size of the Gaussian roughly matches the observed off-axis 
broadening of a point source image. This is not unexpected 
since bright sources are easily detected in any method and, 
for faint sources, Poisson noise removes the importance of 
the exact PSF shape. Assuming a circular Gaussian PSF 
gives higher weight to sources with a central concentration 
of events. This is superior to a cell detect method which 
only looks for an excess of counts in an arbitrary source 
region and does not depend on the distribution of events 
within the region. 

The source detection process was repeated using the 
CXC source detection tool wavdetect (Freeman et al. 
2002). The wavedetect tool is similar to our method in 
that it is more likely to detect sources with a central con- 
densation, wavdetect was applied on spatial scales from 1 
to 16 pixels in logarithmic steps using a significance thresh- 
old for source detection corresponding to less than a 10 -6 
chance probability of detection due to local background 
fluctuations. The results were consistent with our method 
at the equivalent significance level. 

Point-source counts and spectra were extracted from 
within the 95% encircled-energy aperture of the model 
PSF. Background regions were typically chosen from an- 
nular regions surrounding the source regions except in 
crowded regions of the field where background regions 
adjacent to the source were used. The background- 
subtracted counts within the source regions were scaled 
to obtain the aperture-corrected count values. The 
background-subtracted point source detection limit is 12 
counts for the 2.8 minimum S/N threshold and a mini- 
mum 5cr above background. The resulting source posi- 
tions, count rates, and S/N ratios are listed in columns 
2-5, respectively, of Table 2. 

3.2. Counterparts & Source Identifications 

Optical properties of the X-ray sources were deter- 
mined using archival Hubble WFPC2 images, the cat- 
alogue (complete to V < 21) of Perclmuter & Racine 
(1995), and Digitized Sky Survey (DSS) images. Addi- 
tionally, compilations of supernova remanant (SNR) can- 
didates (Matonick & Fesen 1997), HII regions (Petit et al. 
1988; Kaufman et al. 1987), globular clusters (Perelmuter, 
Brodic & Huchra 1995; Chandar et al. 2001), novae (Shara, 
Sandage, & Zurek 1999), and stellar objects (Zickgraf 
& Humphreys 1991; Ivanova 1992; Zickgraf, Szeifert, & 
Humphreys 1996; Sholukhova et al. 1998) were queried for 
spatial correlations to the X-ray sources. Objects within 

6 available from pttp: / /cxc. harvard.edu/caldb/download. html 



the 3<7 uncertainty of the X-ray source positions are con- 
sidered potential counterparts to the X-ray sources. 

Archival Hubble WFPC2 images of portions of the 
Chandra field were searched for potential optical coun- 
terparts to the X-ray sources based on spatial coincidence. 
For most of the images there were V-type filters (F555W, 
F606W, and F547M) available. For a subset of these ob- 
servations there were images with filters corresponding ap- 
proximately to standard U (filter F336W), B (F439W), 
V (F555W), R (F675W), and / (F814W) for each field. 
When there were multiple images, using the same filter, 
they were combined using the IRAF task crrej to remove 
cosmic ray contamination. For a few images in which there 
was a position offset the task xregister was used to cre- 
ate matching images and the task gcombine was used to 
combine the images and to remove cosmic rays. A list 
of point sources detected in each field in a V filter image 
was generated using the star-finding algorithm daof ind in 
DAOPHOT (Stetson 1987). Photometry was performed for 
each source using a 0."3 radius aperture and concentric 
background annuli of inner radius l."0 and 0."5 width. 
For extended sources, a radius as large as l."0 was needed 
for the source aperture. Where multiple optical objects 
fall within the 3tr width of the Chandra model PSF at the 
location of an X-ray source, we selected the most proba- 
ble candidate optical source based on (in order of prefer- 
ence) optical brightness, colors, or distance from the X-ray 
source centroid. 

A total of 66 X-ray sources lie within one or more of the 
Hubble imaging fields. Of these sources, 34 have potential 
optical counterparts based on spatial coincidence with an 
average separation between X-ray source and optical can- 
didate of ^1.5cr. Forty-three of the 66 sources lie in the 
bulge of M81 and 13 of these have optical candidates. In 
contrast, 21 disk sources are within the Hubble fields and 
19 of these have optical candidates based on spatial coin- 
cidence. The search for optical candidates was repeated 
using an artificial distribution of X-ray source positions 
and position uncertainties (while preserving the radial dis- 
tribution of sources). Of 89 artifical X-ray sources, 51 fell 
within one or more of the Hubble imaging fields and 16 
of these had optical candidates (5 candidates in the bulge 
of 36 possible and 11 candidates in the disk of 15 possi- 
ble). Thus, we expect about 20 ±4 of the 34 X-ray sources 
with potential optical counterparts are simply chance co- 
incidence with a greater probability of chance coincidence 
for sources in the disk. Nevertheless, the Hubble observa- 
tions place valid upper limits to the optical luminosity of 
the X-ray sources located within their fields. 

The apparent visual magnitudes of the optical candi- 
dates (from F555W, F606W, or F547M filter measure- 
ments) are listed in column 6 of Table 2. The values have 
not been corrected for charge transfer efficiency (CTE) ef- 
fects nor color-corrected to obtain true Johnson V magni- 
tudes. CTE effects are estimated to be of order 0.01 mag- 
nitude while the color-correction depends on the intrinsic 
properties of the source and on interstellar reddening. The 
correction is typically less than 0.3 magnitudes. For two 
of the sources, numbers 52 and 157, only F814W (/-band) 
images were available and they are denoted as such in col- 
umn 6 of Table 2. Potential optical counterparts to 7 X- 



5 



ray sources appear extended in the Hubble images. Three 
of these are listed as globular clusters in the catalogue of 
Chandar et al. (2001, see below). Two others, source num- 
bers 82 and 127, are probably also globular clusters (see 
below). The remaining 2 sources that appear extended 
in the Hubble images are denoted "Extended?" in column 
10 of Table 2. X-ray sources located within Hubble fields 
but without optical candidates are designated as ">HST" 
in column 6 of Table 2 to indicate they are fainter than 
the Hubble limiting magnitude of ^27 magnitudes. Vi- 
sual magnitudes of objects in the Perelmuter & Racine 
(1995) catalogue within 3er of X-ray sources are listed as 
upper limits in column 6 of Table 2 unless superseded by 
ffubMe-derived values. 

Archived radio observations were used to search for po- 
tential radio counterparts to the X-ray sources. Radio 
data obtained using the Very Large Array (VLA) of the 
National Radio Astronomy Observatory (NRAO) 7 as part 
of an on-going study of SN 1993J were kindly provided 
by N. Bartcl and M. Beitcnholz for our use. For the 
present work, we analyzed a radio map made from a 1999 
November 23 observation at a wavelength of 6 cm (4985 
MHz) in the B configuration. The angular resolution is 
approximately l."5, and the root-mean-square noise level 
is 25 fj,Jy. The image was prepared using standard radio 
data reduction techniques using the NRAO software pack- 
age AIPS (Astronomical Image Processing System). The 
field of view of the radio map is approximately 8/5x8/5 
centered on the optical nucleus of the galaxy. Of the 177 
Chandra-detected X-ray sources, 91 fell within this field 
of view. The two most luminous radio sources in the field 
of view are the nucleus and SN 1993J, and both of these 
sources are known to be time- variable in their radio emis- 
sion. Only two other X-ray sources were clearly detected 
above the 3cr level in this radio data, source numbers 79 
(55±17 /iJy)and 155 (95±26 /iJy). They are designated 
"Radio" in column 10 of Table 2. A source at the location 
of Einstein source X6 was detected at 3.6 cm, as discussed 
Further analysis of the radio properties of the 



4.2 



X-ray sources will be the subject of a separate paper (T. 
Pannuti et al., in preparation). 

There are 5 SNRs in the tabulation of Matonick & Fesen 
(1997) coincident with X-ray sources including Einstein 
source X6. As discussed in §4.2, the X-ray emission from 



X6 is not from a SNR. All four remaining SNR candidates 
are located along spiral arms, as expected for core-collapse 
supernovae from young massive stars, and are listed with 
the designation "SNR" in column 10, Table 2. 

There are no X-ray sources coincident with any of the 
25 globular clusters identified by Perelmuter et al. (1995). 
The positions of four globular clusters in the tabulation 
by Chandar et al. (2001) are within the 3<r error circle of 
X-ray sources as listed in column 10, Table 2, with the 
designation "GC". Three of these, corresponding to X- 
ray source numbers 141, 146, and 158, are located near 
a prominent spiral arm yet the colors reported by Chan- 
dar et al. (2001) place the optical candidates among the 
typical population of old (^1 Gyr) clusters. The remain- 
ing candidate, corresponding to number 148 of Table 2, 
is also an old globular cluster but its projected position 
is between spiral arms. In addition, one of the bright ob- 



jects in the catalogue of Perelmuter & Racine (1995) that 
is coincident with X-ray sources has an optical magnitude 
18<F<21 and colors in the ranges 0.5<(5 - V)<1.1 and 
0.3<(^ — i?)<0.7. This object, source number 127, meets 
the criteria given by Perelmuter et al. (1995) for globular 
cluster candidates and is designated "GC?" in column 10 
of Table 2. Also listed with the same designation are X- 
ray source numbers 67 and 82. Although the optical can- 
didates of both objects exceed the Perelmuter & Racine 
(1995) brightness criteria, Ghosh et al. (2001) have ar- 
gued that they are very likely globular clusters (see also 
$J- 

Three hundred and ninety of the 492 HII regions tab- 
ulated by Petit, Sivan, & Karachentsev (1988) are within 
the ACIS imaging field of view. Twelve of these are coin- 
cident with Chandra x-ray sources. All are located along 
M81 spiral arms. They are designated "HII" in column 10 
of Table 2. 

Three Chandra sources are coincident with foreground 
stars. These are labeled with a designation in column 
10 of Table 2. Two of these, X-ray source numbers 56 
and 60, are the PPM catalogued stars used to determine 
absolute X-ray source positions (Table 1). Inspection of 
DSS images to search for uncatalogued bright star-like ob- 
jects revealed no additional foreground star candidates. X- 
ray source number 8, however, is very likely a background 
galaxy based on its shape in DSS images. 

3.3. X-ray Spectral & Timing Analysis 

Statistically-constrained model fits could be achieved for 
a total of 39 sources in the M81 field. Detailed X-ray prop- 
erties of SN 1993J, Einstein source X6, and the nucleus are 
presented in §0 and of the three brightest supersoft sources 
in Swartz et ah (2002). Best-fit spectral model parameters 
for the remaining 33 bright sources are listed in Table 3. 

The X-ray spectra of these 33 sources were fit to ab- 
sorbed power law models characterized by the photon in- 
dex T; to absorbed Raymond- Smith spectral models rep- 
resenting emission from low density, optically-thin plasma 
characterized by the temperature kT; and to absorbed 
blackbody models. The metal abundances in the ther- 
mal plasma models were assumed to be 3% of their solar 
value consistent with the results of Kong et al. (2000) for 
the galaxy as a whole. Spectra were grouped to contain a 
minimum of 20 counts per spectral bin and fit to models 
using the XSPEC spectral fitting package (Arnaud 1996). 

The power law model provides a significantly better fit 
to all but two of the X-ray sources. The two exceptions, 
X-ray source numbers 160 and 161, are best fit by black- 
body spectral models with effective temperatures 0.2 and 
0.57 keV, respectively. The best- fit absorption column 
density, power-law index (or blackbody temperature), fit 
statistic, and unabsorbed 0.3 - 8.0 keV luminosities are 
listed in columns 2 - 5 of Table 3. Quoted errors are 90% 
confidence limits for a single interesting parameter based 
on the x 2 n t statistic. Absorption column densities were 
constrained in the fitting procedure to be at least as large 
as the Galactic column density along the line of sight to 
M81. This resulted in a best-fit value of Njj equal to this 
lower limit for 8 sources. Therefore, the column densities 
are listed as Njj — 4.0 ± 0.0 in column 2 of Table 3 for 



7 The NRAO is a facility of the National Science Foundation operated under a cooperative agreement by Associated Universities, Inc. 
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these sources. 

The luminosities for the remaining sources in Table 2 
were estimated assuming an absorbed power law spectrum 
with photon index T = 1.5 and hydrogen column density 
jV 2 ci = Nh/10 20 = 11.0 cm~ 2 . These are the average val- 
ues for the 31 sources in Table 3 whose spectra are best 
fit with a power law model. For comparison, the Galac- 
tic hydrogen column density in the direction of M81 is 
N 2Q = 4.0 cm -2 (Stark et al. 1992). The source detec- 
tion limit of 12 counts corresponds to an observed flux of 
-Fo.3-8fceV = 1-9 x 10~ 15 ergs cm~ 2 s" 1 for sources on 
device S3 (F , 3 ^ skeV = 2.5 x 10~ 15 ergs cm~ 2 s _1 on 
FI devices) or an unabsorbed 0.3 - 8 keV luminosity of 
L x = 3.4(4.5) x 10 36 ergs s _1 for sources on the BI (FI) 
devices. The luminosities of all sources are listed in col- 
umn 9 of Table 2. 

The Kolmogorov-Smirnov statistic was used to test 
these sources (including background) for constant count 
rates over the duration of the primary observation. The 
Kolmogorov-Smirnov statistic results are listed in col- 
umn 6 of Table 3. Sources with a low value of Pks have 
a high probability of being variable. 

The 0.3 - 8.0 keV light curves were used to estimate the 
overall X-ray power density spectrum (PDS) of the X-ray 
sources listed in Table 3. Light curves were sampled at the 
3.24104 s nominal ACIS frame time and used to compute 
the Leahy-normalized power spectra of each source. The 
average power in each normalizcd-PDS was 2.0, which is 
consistent with what is expected from Poison noise (Leahy 
et al. 1983). Fluctuations up to 10 - 15 are commonly seen 
in the normalized power spectra with a maximum power in 
one or two frequency bins typically between 15 and 20 for 
all the sources with one exception. Fluctuations of this or- 
der are due to noise. For the one exception, X-ray source 
number 69, the average normalized-PDS power was 2.0 
and the power of all the frequencies was less than 15 ex- 
cept at 0.0245 Hz (corresponding to a 40.8 s period), which 
has a power of 28. Given that there are 4096 frequency 
bins in the PDS for each of 33 sources, the probability of 
seeing one peak with power 28 is 12%. Therefore we con- 
clude that the detection of the peak with power 28 could 
be a statistical fluctuation and that no X-ray pulsations 
are detected in any of the sources with the existing data. 

4. INDIVIDUAL M81 SOURCES 

The eleven brightest sources in the Chandra field ap- 
pear to exceed the Eddington limit luminosity for a 1.5 
Mq spherically-accreting object. It should be noted that 
many accreting compact sources in our Galaxy and in the 
Magellenic Clouds that exceed the Eddington limit are in 
fact neutron stars with episodes of high X-ray luminos- 
ity (see, for example, the compilation in Grimm, Gilfanov, 
& Sunyaev 2001). The four brightest sources in M81 are 
SN 1993J, Einstein source X6, the nucleus, and the bright- 
est supersoft source candidate. The supersoft source is 
discussed in detail by Swartz et al. (2002) and no further 
analysis is given here. Details of the other three sources 
are given in the following subsections follo wed by a brief 
discussion of the remaining bright sources (§4.4). 



4.1. SN 1993 J 

Supernova (SN) 1993J was discovered 1993 March 28 
and was observed by Chandra 2594 days after outburst. 



The pre-supernova star was likely a ~17 M Q star that 
lost all but ~0.2-0.4 M Q of its hydrogen envelope prior to 
explosion (see Wheeler & Filippenko 1996 for an early re- 
view) . This mass was lost through a combination of stellar 
winds and mass transfer to a binary companion, possibly 
involving a common envelope phase (Podsiadlowski et al. 
1993; Woosley et al. 1994; Nomoto, Iwamoto, & Suzuki 
1995). A circumstellar medium (CSM), a relic of this mass 
loss, surrounded the system at the time of explosion as ev- 
idenced by its early radio (Van Dyk et al. 1994) and X-ray 
signatures (e.g., Zimmerman et al. 1994). 

Monitoring of this emission has continued since discov- 
ery. X-ray observations by ROSAT (Immler, Ashenbach, 
& Wang 2001) and ASCA (Uno et al. 2002) began 6 and 8 
days after the explosion and continued intermittently for 
1817 and 564 days, respectively. Radio monitoring con- 
tinues (e.g., Marcaide et al. 1997; Bietenholz, Bartel, & 
Rupen 2001). The evolution of optical line profiles also 
shows an increasing contribution from SN debris interact- 
ing with surrounding gas (Patat, Chugai, & Massali 1995; 
Houck & Fransson 1996; Matheson, et al. 2000). 

The early X-ray and radio data has been successfully 
explained in terms of the circumstellar interaction model 
(Fransson, Lundqvist, & Chevalier 1996). In this model, 
the radial flow in free expansion of the SN debris steep- 
ens into a shock wave in the CSM, heating this gas to 
T ~ 10 9 K. This outer shock cannot be freely expanding 
so the supernova ejecta interacts with the hot gas and a 
second (inner, or reverse) shock front develops. The lower 
temperature (T ^ 10 8 K) reverse shock gas dominates the 
soft xX-ray emission with a flux dependent on the density 
gradient of the SN ejecta and on absorption in the cooling, 
post-shock, gas. 

4.1.1. The Chandra X-ray spectrum of SN 19933 

The observed X-ray spectrum of SN 1993J, the best- 
fit model spectrum (x 2 = 94.6 for 90 dof), and the fit 
residuals are shown in Figure [j]. Fits were applied to the 
spectrum in the energy range 0.3 - 5.0 keV because of 
a lack of signal at higher photon energies, including in 
the Fe K band at ~ 6.5 keV. The model consists of two 
low-temperature absorbed thermal emission-line (vmekal) 
components at 0.35±0.06 keV and 1.01±0.05 keV both 
with N20 = 40.5 ±0.9 cm -2 and a high-temperature mekal 
component with kT = 6.0±0.9 keV, A 20 = 4.9±0.1 cm" 2 . 
Abundances in the low-temperature components are con- 
sistent with solar values with the exception of N (15 times 
solar, though see below), Mg (0.2), Si (1.4), Fe (1.6) and 
He, C, O, and Ne which are all consistent with an abun- 
dance of zero. The high-temperature component is consis- 
tent with subsolar abundances but is not sensitive to this 
parameter. 

The two low-temperture components are needed to 
fit the two peaks present in the observed spectrum (at 
0.78 ± 0.02 and 0.98 ± 0.02 keV, respectively, Figure |). 
These peaks are dominated by Fe L emission from differ- 
ent ionization stages of Fe at the different temperatures. 
In addition, He-like Si and Mg (from the 0.35 keV com- 
ponent) and He- like S, and H-like Si and Mg (from the 
1.01 keV component) are observed in the data. Notably 
absent is the O VIII Lya line at 654 eV. This line is 
predicted to be strong over a range of plasma tempera- 
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tures from -0.1 to 0.5 keV (e.g., Nahar 1999). The high- 
temperature component is needed to account for X-ray 
flux above ~1.5 keV. No emission lines are produced by 
this component. In fact, replacing the mekal component 
with a bremsstrahlung model at the same temperature 
provides an equally-acceptable fit to the spectrum. The 
elemental abundances were allowed to vary in the two low- 
temperature component model fits (with the abundances 
and absorbing columns constrained to be equal between 
the two components). Beginning with solar ratios, the 
abundance of each element was individually varied while 
holding the abundances of the remaining elements fixed 
at their current values until the model converged. This 
procedure was repeated for all a-chain elements from He 
through Ni and for N, Na, and Al. The largest reductions 
in x 2 occurred for N, O, and Fe. A large N abundance 
is required to fit a broad feature near the N VII Lya line 
at 0.5 keV. This feature lies just below the neutral O ab- 
sorption edge at 0.53 keV and may, instead, be an artifact 
of a deeper O edge. This interpretation is also more con- 
sistent with evolutionary models of a massive progenitor 
star (e.g., Thielemann, Nomoto, & Hashimoto 1996) that 
do not predict an over-abundance of N. The model re- 
quires no O consistent with the lack of observed O VIII 
Lya. The over-abundance of Fe needed to fit the spectrum 
may indicate some of the explosively-synthesized material 
has been transported to the outer regions of the SN ejecta 
which are presently entering the reverse-shock region. 
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Fig. 1. — Observed X-ray spectrum of SN 1993J (top) and model 
fit residuals (bottom). Curves trace the full model (solid), 0.35 keV 
(dashed) and 1.01 keV (dot-dashed) vmekal components, and the 
harder, 6.0 keV, mekal component (dotted). 

The combination of low- and high-temperature com- 
ponents is consistent with the standard CSM interaction 
model. The low-temperture emission in this scenario orig- 
inates in the reverse shock region while the hard compo- 
nent comes from the much hotter forward shock. The low- 
temperature components are more heavily absorbed than 
the hard component consistent with an intervening dense 
cool shell of gas at the contact discontinuity between the 
forward and reverse shock fronts. The reverse shock front 
is distorted by Rayleigh- Taylor instabilities (Chevalier & 
Blondin 1995). Emission from regions of differing densi- 
ties (and, hence, cooling rates) will have different char- 
acteristic temperatures. The two low-temperature model 



components may therefore only approximate a range of 
temperatures present in the reverse shock region. 

In contrast, the forward shock is located further from the 
high-density contact region. High-resolution VLBI images 
(e.g., Bietenholz et al. 2001) shows the outer shell to be 
highly circular. While a single temperature forward shock 
may be favored in this case, we note that the temperature 
of the high-temperature component is uncertain due to the 
lack of signal above ~5 keV. The relatively low tempera- 
ture of this component, ^7 x 10 7 K, suggests relatively fiat 
ejecta and CSM density profiles at the time of the Chandra 
observation. 

The average densities in the reverse shock region and in 
the intervening absorbing shell can be estimated from the 
model parameters. The radius of the interaction region 
at 2594 days was approximately 2 x 10 17 cm based on ex- 
trapolating the observed angular size at 1893 days (Bartel 
et al. 2000) to the time of the Chandra observation and 
assuming an average velocity of ^7000 km s _1 (Matheson 
et al. 2000) during the interval. The width of the shell is 
of order 10% of this radius. The number density in the 
thermal emission region is then ^5 x 10 4 cm~ 3 and that 
in the cool absorbing shell is ^2 x 10 5 cm~ 3 . This implies 
the mass in the shell is of order ~0.4 to 2.0 M Q . The lower 
value is consistent with the value of ~0.3 M Q inferred from 
the deceleration observed in radio images of SN 1993J at 
an age of ~5 yr (Bartel et al. 2000). 

Conceivably, an underlying neutron star could also con- 
tribute to the hard emission detected from SN 1993 J. 
While the dynamics of the hot, radioactively-heated, gas 
at the center of the SN is unknown, accretion at a rate 
above a few 10 -8 M Q yr -1 creates substantial X-ray emis- 
sion. However, the emission in the Chandra band will 
be absorbed by the overlying SN debris. Assuming for 
simplicity a uniform-density expanding sphere of gas, the 
ejecta provides a total column density of order 10 22 cm~ 2 
per solar mass of ejecta at the time of the Chandra ob- 
servation. Most of this material is in the form of metals. 
Thus the effective hydrogen column is orders of magnitude 
higher. Woosley, Pinto, & Hartmann (1989) predicted no 
detectable X-ray flux below ~10 keV from an accreting 
neutron star in the center of SN 1987A at an age of 2500 
days with the exception of the Fe Ka fluorescence line. 
A similar conclusion was reached by Xu et al. (1988). 
Though the amount of material ejected by SN 1993J is 
nearly an order-of-magnitude less than in SN 1987A, there 
is as yet no compelling evidence of a neutron star in the 
X-ray spectrum of SN 1993J at an age of ~7 yr. 

4.1.2. The X-ray light curve of SN 19933 

The 0.3 - 8.0 keV luminosity of SN 1993J on day 
2594 was 4.8 x 10 38 ergs s -1 . The flux in the 0.1 - 
2.4 keV ROSAT and 1 - 10 keV ASCA energy bands 
are ~ 3.3 x 10 38 and ~ 3.6 x 10 38 ergs s^ 1 , respectively. 
ROSAT (Immler, Ashenbach, & Wang 2001) and ASCA 
(Uno, et al. 2002) light curves, along with the Chandra 
data, are shown in Figure ||. Visual inspection shows that 
the light curves are not simple power laws but that the 
rate of decline of the X-ray luminosity increases after ~50 
to 100 days. The best-fit broken power law model for the 
ROSAT light curve has its break at ~220 days. The lumi- 
nosity declines as L a £-o.24±o.o4 p r j or ^ tjhe break and 
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as L oc ^- - 62 ± 07 a t later times. The break occurs at 
about 45 days in the ASCA energy range but the change 
in slope is less pronounced, evolving from L oc £-°- 57 ± - 78 

to L oc <-0.84±0.24^ 
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Fig. 2.— Observed X-ray light curve of SN 1993J in the ROSAT 
(top) and ASCA (bottom) energy bands. The best-fit broken power 
law models are shown as solid lines. 

As originally conceived, the self-similar form of the in- 
teraction model of Chevalier (1982; Fransson, Lundqvist, 
& Chevalier 1996) applies only to the early phases of su- 
pernova evolution when both the SN ejecta and CSM den- 
sity profiles can be represented by single power laws in 
radius. Extrapolation to ~7 yr is inappropriate because 
as the reverse shock progresses through the ejecta the den- 
sity profile flattens while the forward shock has traversed 
some 10 4 yr of previous mass-loss history (assuming a wind 
velocity of order 10 km s _1 and a shock velocity of or- 
der 10 km s _1 ). Furthermore, the importance of radia- 
tive losses behind the reverse shock and departures from 
electron-ion equipartition in the forward shock make nu- 
merical calculations necessary. 

Suzuki & Nomoto (1995) have performed hydrodynamic 
calculations based on models of the SN ejecta that ac- 
curately reproduce the observed optical light curve of 
SN 1993J (Nomoto et al. 1995). Suzuki & Nomoto (1995) 
followed the evolution for ~1000 days and find reasonable 
agreement with X-ray observations for the first 50 days or 
more. Beyond that time, Suzuki & Nomoto (1995) require 
a combination of a steeper CSM density gradient and a 
clumpy CSM morphology to sustain a level of X-ray flux 
comparable to that observed. However, their models pre- 
dict rapid increases in X-ray flux, followed by declines on 
timescales of order months, as the forward shock sweeps 
through individual clumps. This is in contrast to the ob- 
served steady decline of the X-ray light curve over the 
entire monitoring sequence. 

While beyond the context and scope of the present work, 
X-ray observations of SN 1993J warrant further investi- 
gation. As pointed out by Immler et al. (2001), X-rays 
from the interaction region trace the pre-SN evolution of 
the progenitor system and X-rays are a direct means of 
accessing this evolution in detail as has been undertaken 
recently for SN 1987A (Park et al. 2002). 



4.2. Einstein Source X-6 

The brightest non-nuclear source in the Chandra field 
is the Einstein-discovered source X6 (Fabbiano 1988) lo- 
cated ~1' to the southeast of and along the same promi- 
nent spiral arm containing SN 1993 J. The X-ray flux in 
the Einstein observation was 9.5 x 10~ 13 ergs cm~ 2 s _1 , 
placing X6 in the class of Ultra-Luminous X-ray sources 
(ULXs) whose luminosities (^10 39 ergs s _1 ) far exceed the 
Eddington limit for spherically-accreting ~1.5 M Q objects. 
The X-ray flux from X6 has remained remarkably steady 
throughout its observed history. X6 is coincident with a 
weak radio source (Fabbiano 1988) and with an optically- 
identified SNR candidate (Matonick & Fesen 1997). Chan- 
dra can resolve sources on scales smaller than the 90 pc 
(5") diameter reported for the SNR. 

4.2.1. The Chandra X-ray spectrum of X6 

The high X-ray flux from X6 leads to a pileup of events 
in the ACIS detector and to statistically-significant de- 
tection of events during frame transfer. Events detected 
during frame transfer appear as a streak or trail through 
the source along the detector readout direction. These 
events are not piled up and therefore represent the true 
count rate (when properly time-scaled) and source spec- 
trum although spread over a large spatial region. The 
0.3 - 8.0 keV count-rate of X6 is 0.21 c s" 1 compared to 
the 0.38 ± 0.04 c s~ x rate deduced from the readout trail. 
The spectrum was fit using two identical models but with 
one model convolved with the pileup model developed by 
Davis (2001) as implemented in XSPEC v. ll.l.Ou. This 
combination reflects contributions from pileup and non- 
pileup spatial regions: Piled-up events are localized to the 
central few pixels containing the majority of the detected 
counts. The spectral extraction region includes this cen- 
tral region and additional source counts from surrounding 
pixels. A linear energy grid is necessary to apply the con- 
volution model. Therefore, the spectral analysis is per- 
formed with the Level 2 event data using CXC-provided 
response and ancillary response files instead of those pro- 
vided by L. Townsley which utilize a piece- wise linear en- 
ergy grid. As a test of our procedures, we repeated the 
analysis of S5 0836+710 undertaken as a demonstration 
by Davis (2001) and derived model parameters consistent 
with the values presented in that paper. 

Spectral models were applied to photons from X6 in the 
0.3 to 10.0 keV range using a 2" extraction region. The 
extension to 10 keV was made because there are substan- 
tial source counts above 8 keV and the observed spectrum 
shows a flattening above ~7 keV characteristic of pileup 
(Figure ^), a feature helpful for constraining the model 
parameters. The 2" extraction region is the same size re- 
gion used by Davis (2001) in his analysis of S5 0836+710. 
The fitting procedure resulted in a slight adjustment of 
the pileup parameters from the values reported in Davis 
(2001). The event grade morphing parameter, a, and the 
fraction of the extraction region not experiencing pileup, 
1 — /, both must be increased slightly because the PSF is 
slightly asymmetric and broadened at the ~1' off-axis posi- 
tion of X6 relative to the on-axis location of S5 0836+710. 
Inspection of the high-resolution image of the model PSF 
shows the asymmetry increases the probability that the 
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second photon in a two-photon event will enter an adja- 
cent pixel rather than a corner pixel and hence will be 
registered as a good grade. The broadening of the PSF 
results in a larger fraction of the detected photons falling 
in the wings of the PSF where pileup does not occur. The 
resulting pileup parameters are a = 0.585 compared to 0.5 
used by Davis (2001) and 1-/ = 0.063 (compared to 0.05) 
as determined from the best-fit spectral model. 
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Fig. 3. — Observed X-ray spectrum of Einstein source X6 (top) 
and model fit residuals (bottom). The curve traces the best-fit disk 
blackbody model spectrum with pileup modeled as in Davis (2001). 

The best-fit model for the observed spectrum of X6 is 
an absorbed disk blackbody. Makishima et al. (2000) have 
modeled the ASCA spectrum of X6 and other ULXs us- 
ing the disk blackbody model they developed (Mitsuda 
et al. 1984) for modeling the high soft state of accreting 
black holes. X-ray emission in this state originates from 
an optically thick accretion disk and the model is basi- 
cally a superposition of blackbody emission from different 
disk annuli with local disk temperatures scaling with the 
disk radius as i?" 3 / 4 . Thus the model spectrum in the 
X-ray regime is dominated by the innermost disk temper- 
ature, Ti n , with a normalization scaling as the disk ge- 
ometry; K oc (R in /D) 2 cos(9) where Ri n is the innermost 
disk radius, D the source distance, and 9 the disk incli- 
nation. This model provides a fit statistic \ 2 = 295.3 
for 283 dof. The best-fit innermost disk temperature 
is T in — 1.03 ± 0.11 keV, the corresponding radius is 
Rin = 161 ± 16 km, and the absorbing column density is 
A^o = 21. 7± 1.0 cm" 2 . The disk blackbody model param- 
eters correspond to an 18 M accreting, non-rotating, ob- 
ject and a bolometric luminosity of ~2.7x 10 39 ergs s _1 ac- 
cording to the relations given by Makishima et al. (2000). 
This luminosity is equivalent to the Eddington luminos- 
ity for the derived mass. The parameters obtained from 
analysis of the ASCA spectra (Mizuno 2000, Makishima 
et al. 2000) are T in = 1.48 ± 0.08 keV, R in = 83 ± 8 km, 
and A^o = 21 ± 3 cm -2 . The higher temperature implies 
a lower mass since Ti n oc M -1 ' 4 (eq. 12, Makishima et 
al. 2000) but a higher luminosity L oc R 2 n Ti nl implying 
the luminosity of X6 exceeds the Eddington limit. Mak- 
ishima et al. (2000) argue that, if the compact object is 
a Kerr black hole, then the inner radius can be reduced 
by as much as a factor of 6 thereby reducing the derived 



bolometric luminosity to values below the Eddington limit 
for the temperature-estimated mass. In contrast, the fit 
parameters to the Chandra data leads to self-consistent 
values of bolometric luminosity, mass of the compact ob- 
ject, and the associated Eddington limit luminosity. 

Mizuno (2000) analyzed 7 individual ASCA observations 
of X6 and found a temperature Ti n — 1.3 ±0.1 keV in two 
of the observations and 1.6 ± 0.1 keV in the remaining ob- 
servations. The value derived here, T m = 1.0 ±0.1 keV, is 
significantly lower than these values. It is unclear if this 
is a real effect or is an artifact of the pileup model. 

A power law model fit to the Chandra spectrum was 
statistically less acceptable. The best-fit result with a and 
1 — / fixed as above is \ 2 — 342.8 for 284 dof; for an ab- 
sorbing column iV 2 o = 39.5±2.0 cm" 2 , and a photon index 
T = 2.1 ± 0.1. Across the 90%-confidence range of a de- 
termined from the disk blackbody model, 0.56<a<l, the 
range of acceptable power law indices is 2.0<T<2.3. The 
largest systematic contributions to \ 2 f° r this model occur 
just above the 2 keV Ir-M edge. This is precisely where 
contributions from pileup from photons at the peak of the 
energy distribution (at ~1.0 - 1.5 keV) occurs. Adding a 
broad Gaussian line to the model significantly improved 
the fit to x 2 = 298.0 for 281 dof. The resulting absorbing 
column is N20 = 31.2 ± 2.0 cm" 2 and the photon index is 
T = 1.71 ± 0.09. The aperture-corrected model predicted 
flux in the 0.3 - 8.0 keV band is 3.9 xlO -12 ergs cm" 2 s" 1 
corresponding to a luminosity of L = 6.0 x 10 39 ergs s _1 . 

An optically-thin thermal plasma model can also pro- 
duce an acceptable fit provided the abundance is kept 
low to reproduce the observed continuum-dominated spec- 
trum. The best-fit parameter values for X6 are kT = 
3.5 ± 0.4 keV, N 20 = 30.7 ± 1.8 cm" 2 , for X 2 = 310.0, 
283 dof, and the metal abundance constrained to 0.03 of 
the solar value. This model produces few strong spectral 
lines with the exception of Fe Ka at 6.7 keV. Adding a 
line at this energy to either the disk blackbody or power 
law models is not statistically significant. 

The results given above are consistent with another test 
we conducted: The spectrum of X6 was extracted from 
an annulus with a l."5 inner radius (compared to the 2" 
outer radius used for applying the pileup model). In this 
case pileup is not an issue and a simple absorbed disk 
blackbody or power law model is sufficient. The result- 
ing model parameters were consistent with those quoted 
above though the formal errors are considerably larger be- 
cause of the lower number of source photons detected in 
the extraction region. 

4.2.2. X6 Radial Profile 

A two-dimensional Gaussian model fit to the spa- 
tial distribution of X-ray events places source X6 at 
9 fl 55 m 32.98±0.08 s , ±69°0'33.4±0.4" (ignoring the abso- 
lute uncertainty in source positions, see ^ and Table 1). 
Source X6 is coincident with a large SNR candidate. Wang 
(1999) identified several X-ray-bright sources apparently 
associated with SNRs based on ROSAT observations of 
M101 and concluded that the blast wave energies of these 
SNRs exceed theoretical predictions of supernova explo- 
sions by factors of 30 or greater. In the case of X6, the 
superb angular resolution of the Chandra image can be 
used to determine the extent of the source on scales much 
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less than the ~5" size of the candidate SNR (No. 22, Ta- 
ble 10, Matonick & Fesen 1997). 

Figure f| displays the radial profile of source X6 along 
with a model of the radial profile for a point source located 
at the off-axis position of X6 and the observed radial pro- 
file of SN 1993J. The total number of counts in the model 
PSF and the SN 1993J profile were scaled to 1.6 times the 
total number of counts detected in the X6 profile. This 
accounts for the 40% pilcup estimated for source X6. As 
can be seen, the majority of the pileup occurs in the in- 
nermost two radial bins (with an area of 9 pixels). The 
profile of SN 1993J is slightly more concentrated than ei- 
ther the model PSF or the X6 profile consistent with its 
closer proximity to the aimpoint (20" compared to ~58" 
for X6). There is no evidence that X6 is an extended 
source; including in the 5" region occupied by the SNR. 
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Fig. 4. — Radial profile of source X6 (solid). X6 is located ~1' 
off-axis. The profiles of SN 1993J (dotted; ~20" off-axis) and the 
1.5 keV, l'-off-axis, model PSF (dashed) are shown for comparison. 
There is no evidence for source extension in the X-ray profile of X6. 



4.2.3. Potential X6 Counterparts 

A weak uncatalogued radio source visible in the 21 cm 
map of Bash & Kaufman (1986) is coincident with the 
Einstein High Resolution Imager position of source X6 ac- 
cording to Fabbiano (1988). The source is not seen in a 
6 cm VLA image taken on 1999 Nov 23 above the ~ 80/iJy 
3cr limit but is detectable in a 3.6 cm image obtained 1994 
Dec 23 at a flux density of ~ 95/xJy, just above the 3cr 
signal-to-noise limit. The radio source extension cannot 
be reliably measured at this low signal level. 

Figure g displays an archival Hubble WFPC2 F555W 
image of the region containing X6 with the X-ray source 
position identified by a l"-radius circle. An optical source 
is clearly present within ~ 0."2 of the X6 location. The 
observed Hubble magnitudes are F336W= 22.8 ± 0.2, 
F439W= 24.1 ± 0.3, F555W= 24.1 ± 0.1, F675W= 23.9 ± 
0.2, and F814W= 23.7 ± 0.5. Estimating the color ex- 
cess from the hydrogen column density obtained from 
the X-ray spectrum, the optical properties are consis- 
tent with an early-type main sequence star of spectral 
class 09 - Bl though the source is relatively bright in 
the i?-band (F675W) image, perhaps due to Ha emis- 



sion. Be stars with circumstellar disks can emit strong 
Ha under certain circumstances. The equivalent Ha flux 
is ~ (5.4 ± 0.9) x 10~ 15 ergs cm~ 2 s _1 assuming all the 
i?-band flux is from this emission line. Matonick & Fesen 
(1997) report the Ha flux from the SNR candidate to be 
1.1 x 10~ 14 ergs cm -2 s -1 . The point source can, there- 
fore, contribute as much as one-half of this amount. Note 
that the point source flux estimate is not contaminated 
by any underlying extended emission because the source 
a nd b ackground extraction regions used in our analysis 
(§3.2) lie wholly within the region occupied by the candi- 
date SNR. 





Fig. 5. — Hubbie/WFPC2 image of the region surrounding source 
X6 taken with the F555W filter. The l"-radius circle at the lower 
left denotes the X-ray position in the Chandra data. The arrow indi- 
cates the direction of north and the line segment denotes east. The 
optical properties of the object at the location of X6 is consistent 
with a main sequence star of spectral class 09 — Bl. 

If the X-ray emission is from an accretion disk then it 
is possible that some or all of the optical emission is also 
from an accretion disk instead of from a companion star 
(or from another object in the field). The color index, 
£ = Bq + 2.51ogFx where Bq is the reddening-corrected B 
magnitude and Fx is the 2-10 keV X-ray flux in /iJy, was 
introduced by van Paradijs & McClintock (1995) to quan- 
tify the observed relationship between X-ray and optical 
flux from accretion-powered X-ray binaries. The average 
value of £ for systems in which the secondary star does not 
contribute significantly to the optical brightness (namely, 
those with low-mass donor stars) is 21.8 ± 1.0 (errors are 
1 standard deviation) . This can be compared to X6 where 
£ ~ 22.5 ± 1.5 (where the error includes an uncertainty 
of 0.15 magnitudes in the color excess). The average col- 
ors for these systems are (B — V) = —0.09 ± 0.14 and 
(U-B) = -0.97±0. 17 compared to (B-V) = -0.1±0.3 
and (U-B) = -1.4 ± 0.3 for X6. 
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The observed X-ray-to-optical flux ratio, ^900, is far too 
large for typical foreground objects. Normal stars have ra- 
tios in the range 1CP 4 to 0.1 (Maccacaro et al. 1988) and 
cataclysmic variables (CVs) have ratios in the range 0.1 to 
10.0 (Bradt & McClintock 1983). The only extragalactic 
sources with very high ratios of X-ray to optical flux are 
the BL Lac sources with ratios in the range 10 to 50 (Mac- 
cacaro et al. 1988). To produce the ratio observed for X6 
requires an optical extinction greater than 4 magnitudes, a 
possibility excluded by the X-ray-measured hydrogen col- 
umn density of ~ 2 x 10 21 cm~ 2 . 

4.2.4. X6 Variability 

Tests for source flux variability were applied to counts 
extracted from both the entire source region and to counts 
in the inner 3x3 pixel region and to the outer source region. 
There is no evidence of pulsations or other variability in 
any of these regions. However, periods of higher flux would 
incur more pileup. This tends to smooth the light curve 
by reducing the observed count rate during these high flux 
periods. 

We performed similar analysis of the extensive set of 
ROSAT observations and found no large-scale variability 
on timescales as short as the ROSAT orbital period. Imm- 
ler & Wang (2001) report a factor-of-two change in the X6 
ROSAT PSPC count rate over a 6 day period in 1993 Nov. 
Our analysis of this data find a modest, maximum 22±7%, 
change in count rate but the data are statistically con- 
sistent with no variation when compared to SN 1993J. 
Mizuno (2000) reports that variability of X6 in the ASCA 
data cannot be assessed because of the unavoidable con- 
tribution from the bright nucleus. 

4.3. The Galactic Nucleus 

The nucleus of M81 has long been an object of study. It 
is optically classified as a low-ionization nuclear emission 
line region (LINER) (Ho, Filippenko, & Sargent 1996). 
The dominant source of energy in LINERs may be me- 
chanical heating by shocks, photoionization by hot stars, 
or photoionization by a low-luminosity active galactic nu- 
cleus (AGN). The nucleus of M81 shows evidence of a low- 
luminosity AGN including a compact radio core (Bieten- 
holz et al. 2000), broad Ha emission (Peimbert & Torres- 
Peimbert 1981), a UV-bright continuum (Devereux, Ford 
& Jacoby 1997) with broad, AGN-like emission lines (Maoz 
et al. 1998), and a power law X-ray continuum (Ishisaki 
et al. 1996, Pelligrini et al. 2000). These are all consis- 
tent with the presence of a 6 x 10 7 M Q object, as inferred 
from dynamical studies (Bower et al. 2000), at the galactic 
center. 

There have been numerous X-ray studies of the nucleus 
including Einstein (Elvis & van Speybroeck 1982), GINGA 
(Ohashi et al. 1992), BBXRT (Petre et al. 1993), ROSAT 
(Radecke 1997), ASCA (Ishisaki et al. 1996, Iyomoto & 
Makishima 2001), BeppoSAX (Pellegrini et al. 2000), and 
XMM-Newton (Page et al. 2002). A summary of the cur- 
rent and previous X-ray observations of the M81 nucleus 
is given in Table 4. In addition to a power law with a 
slope of 1.9, similar to those of luminous Seyfert 1 nuclei 
(Turner & Pounds 1989, Nandra et al. 1997, Terashima et 
al. 2002), low-resolution X-ray observations suggest a soft 
thermal component is present in the nucleus. While both 



advection-dominated accretion flows (ADAFs) and accre- 
tion disk coronae can provide the necessary Comptonizing 
medium to produce the observed power law X-ray spec- 
trum, the source of the thermal component remains an 
outstanding issue. This is also true in other wavebands. 
The UV continuum is weak relative to X-rays and the "big 
blue bump" is absent in M81; perhaps a manifestation of 
a low accretion rate (Ho et al. 1996) or the presence of an 
ADAF (Quataert et al. 1999). 

Table 4 

X-Ray History of the M81 Nucleus 



Date 


Obs. 


■pa 


K 


F c 


Ref 


1979 Apr 


Einstein 




68±58 


~0.1 


1 


1985 Feb 


EXOSAT 


2.1±0.2 


19±9 


2.65±0.25 


2 


1987 May 


GINGA 


2.2±0.2 


60±30 


2.4±0.3 


3 


1990 Dec 


BBXRT 


9 9+O.3 


41±10 


3.6±0.3 


4 


1991 - 94 


ROSAT 


2.5±0.3 


7.4±0.4 


1.24 


2 


1993 - 99 


ASCA 


1.85±0.04 


10.0 


3.50±0.12 


5 


1998 Jun 


SAX 


1.86±0.03 


12 n+ 3 ' 6 

^—3.0 


3.8 


6 


2000 May 


Chandra 


1.98±0.08 


9.4±2.0 


3.25±0.15 


2 


2001 Apr 


XMM 


1.94±0.06 


3.4±0.8 


0.93±0.06 


7 



"Power law photon index 

b Solar-abundance absorption column, in units of 10 20 cm 

c 2 - 10 keV observed flux, in units of 10 -11 ergs cm~ 2 s - 

References.- (1) Fabbiano 1988 (2) this work 

(3) Ohashi et al. 1992 (4) Petre et al. 1993 

(5) Iyomoto & Maikshima 2001 (6) Pellegrini et al. 2000 

(7) Page et al. 2002 



Though past X-ray observations have been unable to 
resolve the nucleus from the surrounding diffuse emission 
and point-like sources, Tennant et al. (2001) have shown 
that some ^10 39 ergs s _1 , or a few percent, of the X-ray 
luminosity in the nuclear region actually originates within 
a rather extended (~ 2/5) region of the bulge. A similar 
conclusion has been reached more recently by Immler & 
Wang (2001) from analysis of ROSAT data and by Page 
et al. (2002) from XMM-Newton RGS spectral analysis. 

The high X-ray flux from the nucleus leads to a severe 
pileup of events in the ACIS detector, making the point- 
source spectrum unsuitable for study (although J. Davis 
has successfully analyzed the nuclear spectrum; J. Davis, 
2002, private communication). In this work, the nuclear 
spectrum was instead extracted from two 10"-wide rectan- 
gular regions spanning the readout trail and offset >30" 
from the nucleus to avoid contamination by events in the 
wings of the nuclear PSF. The background was extracted 
from 4 similar rectangular regions adjacent to the readout 
trail. 

A simple absorbed power law model provides a statisti- 
cally acceptable fit to the nuclear spectrum (\ 2 = 114 for 
104 dof; Figure ||). The resulting power law fit parameters 
are listed in Table 4. The addition of a thermal compo- 
nent did not improve the fit significantly (A% 2 = 0.4 for 
one additional parameter with the temperature parame- 
ter constrained to kT — 0.5 keV). Detection of Fe Ka 
emission could not be confirmed because of the lack of 
source counts above ~5 keV in the extracted spectrum. 
This line was previously detected in Ginga (Ohashi et al. 
1992), ASCA (Ishisaki et al.1996), BeppoSAX (Pelligrini 
et al. 2000), and, recently, XMM-Newton PN (R. Soria, 
private communication) observations of M81. 
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Fig. 6. — Spectrum of the M81 nucleus extracted from the readout 
trail. Shown are the contributions from a T = 1.98 power law (solid 
line) and a weak thermal component (dot-dashed line). The addi- 
tion of the thermal component does not improve the fit significantly 
(A% 2 = 0.4 for 1 additional parameter). 

The 90% upper limit to the thermal model normaliza- 
tion corresponds to a maximum thermal contribution to 
the nuclear flux of 2.8% or a luminosity in the 0.3 - 8.0 keV 
band of ~9 x 10 38 ergs s~ 1 . While this is only a small frac- 
tion of the nuclear emission it is comparable to the total 
luminosity from unresolved emission in the bulge . The 
lack of a significant thermal component in the Chandra nu- 
clear spectrum is in contrast to results based on low spatial 
resolution observations including ASCA, BeppoSAX, and 
XMM-Newton that find one or more thermal components 
improve the fit significantly. In fact, it is argued in ^ that 
any thermal emission from the nuclear region is consistent 
with an extrapolation of the emission from the region be- 
yond ~10" from the nucleus and therefore is unrelated to 
the nucleus. 

4.4. Other X-ray-bright M81 Sources 

In addition to SN 1993J, X6, the nucleus, and the bright- 
est supersoft source (Swartz et al. 2002), there are seven 
X-ray point sources with luminosities exceeding the Ed- 
dington limit for a 1.5 M Q accretor in M81. Four of these 
are the Einstein sources X7, X4, X10, and X3 in order of 
decreasing brightness. These have all also been detected 
in the ROSAT data (Table 2). In addition, source num- 
ber 83 was detected by ROSAT. The two remaining bright 
sources are confused with nearby sources in the ROSAT 
and Einstein observations: Source number 116 is only 16" 
from the nucleus and source number 57 is 19" from (and of 
comparable luminosity to) X10. Thus, aside from the two 
confused sources, the brightest sources must be relatively 
steady sources. 

The X-ray spectra of these sources are unremarkable 
with the exception of source X4 (number 86) which is ex- 
tremely flat, r = 1.08 ± 0.13, and source number 57 which 
has a relatively high column density, N20 = 99.1 ±40. Our 
analysis of the ROSAT PSPC spectra of these two sources 
is consistent with these values (though source 57 is con- 
fused with source X10, number 52). All seven sources are 
within the D25 isophote of M81. Three are within the 
bulge as is the brightest supersoft source (number 132). 



There are no discernable optical or UV counterparts to 
any of these seven sources with the exception of source 
number 146, Einstein source X7, which is coincident with 
globular cluster number 63 of Chandar et al. (2001). 

While the majority of these bright sources are persis- 
tent sources, two Chandra sources of sub-Eddington lu- 
minosity are known to have been much more luminous in 
the past. These are Chandra source numbers 67 (Ein- 
stein source X2) and 82. Ghosh et al. (2001) analyzed 
source number 82 in detail and noted that both it and 
X2 are spatially-coincident with optically-bright objects 
which may be bright M81 globular clusters. The bright 
transient source number 82 reached a peak unabsorbed 
luminosity of ~ 7 x 10 38 ergs s _1 during ROSAT observa- 
tions but was observed at 1.7 x 10 37 ergs s -1 in the Chan- 
dra data (Ghosh et al. 2001). Source X2 is not known to 
have exceeded ~ 2 x 10 38 ergs s -1 . All ROSAT sources 
in the tabulation of Immler & Wang (2001) falling within 
the Chandra field of view have Chandra-detected counter- 
parts (Table 2). No other bright transients are present in 
these two tabulations based on a comparison of the Chan- 
dra count rates derived in this work and the ROSAT count 
rates tabulated by Immler & Wang (2001). 

5. THE M81 BULGE 

The inner Lindblad Resonance, located at an 
inclination-corrected radius of 4 kpc, separates the bulge 
and inner disk from the spiral arms and outer disk of M81 
(e.g., Kaufman et al. 1989, Reichen et al. 1994). This value 
is larger than the 2.5 kpc radius adopted by Tennant et 
al. (2001) in their discussion of bulge and disk emission 
(and based on the observed X-ray morphology). A 4 kpc 
radius circle in the plane of the galaxy corresponds to a 
7/64 x 3/94 ellipse on the plane of the sky with major axis 
oriented at PA 149°. 

The Chandra image shows the bulge X-ray sources con- 
centrated towards the galactic center and excess or unre- 
solved X-ray emission extending away from the nucleus. 
Analysis of this emission is made difficult by the bright 
nucleus. While the high angular resolution of the Chan- 
dra mirrors concentrates most of the nuclear X-rays into 
the central few pixels, the small fraction incident in the 
wings of the PSF accounts for a large percentage of the 
total number of X-rays detected at larger radii. 

In this section, the spatial distribution of the bulge X- 
ray emission is analyzed and compared to the distribution 
observed in other wavebands. The spectral properties are 
then assessed to determine the possible contributions from 
different sources of X-ray emission as a function of position 
within the bulge. 

Of particular interest is the central ~30 - 50" inner re- 
gion or core of the galaxy. Unlike the smooth distribution 
of optical light, there is filamentary Ha emission (Dev- 
ereux, Jacoby & Ciardullo 1995) and excess UV emission 
(Hill et al. 1992, Reichen et al. 1994) in the core. The ori- 
gin of this emission is uncertain. It has been attributed to 
recent star formation activity (Devereux, Jacoby & Cia- 
rdullo 1995) but more likely originates from hot evolved 
post-AGB stars (O'Connell et al. 1992, Devereux, Jacoby 
& Ciardullo 1995). High resolution Hubble images rule 
out massive OB stars as a source of ionization but ioniza- 
tion by shocks originating from nuclear activity remains 
a viable alternative (Devereux, Ford, & Jacoby 1997), a 
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conclusion also reached by Greenwalt et al. (1998). 

5.1. Spatial Distribution 

5.1.1. Broad-band X-ray Surface Brightness 

The unresolved X-ray surface brightness is obtained by 
removing the detected sources and the nuclear readout 
trail from the image. The radial profile of this emission is 
shown in Figure H. The profile asymptotically approaches 
the background level of 0.04 cts pixel -1 at a radius of ~2' 
(~2 kpc). Also shown is the estimated contribution from 
the nuclear point source and the resulting profile with the 
nuclear contribution subtracted. The contribution from 
the nucleus was estimated using a model PSF appropri- 
ate for the observed off-axis location of the nucleus and 
for 1.5 keV photon energies. This energy is above the 
^ 0.8 keV peak of the observed spectrum of the nucleus 
(§4.3) but is representative and avoids interpolation be- 
tween PSF models at different discrete energies. The ra- 
dial profiles show a strong nuclear contribution within the 
central ~30" but which falls much more steeply than the 
observed excess at larger radii. Toward the center of the 
galaxy, the estimated contribution from the nucleus ex- 
ceeds the observed surface brightness due to pileup of nu- 
clear photons. Analysis is therefore confined to the region 
beyond ~10" from the nucleus. 
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Radius (arcsec) 

Fig. 7. — Curves tracing the radial profile of the observed X-ray 
surface brightness (dotted), the model point spread function of the 
nucleus (dashed), and the difference of the two (solid). Errors have 
been omitted for clarity. The background level is ~ 0.04 c s _1 

The excess X-ray emission (with the estimated nuclear 
PSF contribution subtracted) was examined to determine 
the shape of the emission and to compare to profiles at 
other wavelengths. Both the bulge and the underlying 
galactic disk can contribute to the excess X-ray surface 
brightness. If the excess is confined to the disk, then 
its distribution should appear elongated in an elliptical 
pattern consistent with the known inclination of M81. If 
the emission is from the bulge, then its distribution on 
the plane of the sky should appear azimuthally symmetric 
about the nucleus. 

Azimuthal profiles extracted over a range of spatial 
scales show no significant departures from azimuthal sym- 



metry with one exception: There is a slight (~30%) in- 
crease in X-ray surface brightness in a region located 30" 
to 45" to the northeast of the nucleus. This is also the lo- 
cation of a nonthermal highly polarized radio arc that may 
be a small-scale nuclear radio lobe (Kaufman et al. 1996). 
A two-dimensional Gaussian model fit to the excess X- 
ray surface brightness shows a slight elongation along PA 
153±2.4° and an eccentricity of 0.56. This is much less 
than expected from the galaxy disk inclination (eccentric- 
ity 0.86) and slightly smaller than the ~0.66 measured in 
the inner 0.1 - 1 kpc from UBVR isophotes (Tenjes, Haud, 
& Einasto 1998). Qualitatively, this suggests the excess X- 
ray surface brightness is emitted predominantly from the 
spherical bulge but that the disk also contributes. 

A more quantitative estimate of the disk and bulge con- 
tributions can be made by fitting the radial profile of the 
excess X-ray surface brightness with a generalized expo- 
nential of the form S(r) = Er ) e~' r/ '' l ) 1/ " (e.g., de Jong 
1996). The PSF-corrected and background-subtracted X- 
ray surface brightness profile is shown in Figure || along 
with the best fit curves for n = 1, 2, and 4. None of 
the fits are statistically acceptable, with x 2 values ranging 
from 160.4 to 203.1 for 116 dof. Allowing n to vary im- 
proved the x 2 value only slightly to 154.6 for the best-fit 
value n = 1.56 ± 0.30. This value is intermediate between 
an exponential disk (n = 1) and a de Vaucouleurs bulge 
(n = 4 or i? 1 / 4 law) profile. Fitting only the region beyond 
1' of the nucleus reduced x 2 substantially for n = 1 (65.4 
for 65 dof) and results in a best-fit value of n = 1.1 ± 0.1 
(x 2 = 65.3 for 64 dof). The X-ray emission at large radii 
therefore generally follows the exponential disk profile typ- 
ically seen in optical light from spiral galaxies but steepens 
to a bulge-dominated profile near the nucleus. 
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Fig. 8. — Background-subtracted and PSF-corrected excess X-ray 
surface brightness radial profile. Curves represent best fits using a 

generalized exponential function of the form S(r) = Erje - ( r /*0 ' 
with n = 1, 2, and 4. (Curvature increases with increasing n; the 
exponential, n = 1, is the straight line). 
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Fig. 9. — Background-subtracted and PSF-corrected excess X-ray 
surface brightness radial profile (triangles) with optical (squares), 
near-UV (circles), and far-UV (crosses) profiles. Non-X-ray profiles 
have been scaled vertically to match the X-ray data at a radius of 
10" for comparison. 



Fig. 10. — Background-subtracted and PSF-corrected excess X- 
ray surface brightness radial profile (triangles) with histogram dis- 
tribution of X-ray counts detected in resolved X-ray sources (in units 
of c kpc -2 ) overplotted. Profile of resolved source counts has been 
scaled vertically by 2.5 X 10 -4 . The high bin at 50" is due mostly 
to a single source contributing ~4000 X-ray counts. 



Figure |^ compares the X-ray surface brightness profile 
to those observed at other wavelengths. The X-ray surface 
brightness follows the optical (courtsey P. Tenjes 2002, pri- 
vate communication) and near- UV 8 (centered at 2490 A) 
light from old bulge stars out to ~80" but declines more 
rapidly at larger radii. In contrast, the far-UV (1520 A) 
profile falls rapidly to the background level beyond 50" 
of the nucleus (Hill et al. 1992). The resulting UV color 
gradient has been interpreted as a gradient in metal abun- 
dance in the inner ~50" (O'Connell et al. 1992). There is 
no indication of a similar enhancement in the X-ray data 
but, again, the morphology within the inner ~10" cannot 
be determined owing to the X-ray brilliance of the nucleus. 



5.1.2. Resolved Sources 

The radial distribution of the observed counts (exclud- 
ing the nucleus) in the resolved X-ray population super- 
posed on the excess X-ray surface brightness profile is dis- 
played in Figure [l0|. The profiles are remarkably similar 
with the exception of the single high bin due to a single 
bright (supersoft) source (number 132) located ~50" from 
the nucleus. 

The radial distribution of the 53 resolved bulge X-ray 
sources (per unit area) is flatter than the surface brightness 
profile. There is, however, a bias in the source detection 
efficiency in the center of the galaxy where the excess X- 
ray surface brightness is highest and acts as an increased 
background for source detection. This naturally causes a 
flattening in the resolved source number distribution. 



5.2. X-ray Spectrum of the Unresolved Component 

The unresolved X-ray emission may be composed of un- 
resolved point sources and of diffuse, shock-heated, gas. 
The spectrum of the unresolved X-ray emission was ex- 
tracted from 15"-wide annuli centered on the nucleus with 
the inner radius of the first annulus equal to 10". At least 
three spectral components are included in spectral fits for 
each annulus: a nuclear contribution, a power law, and a 
thermal component. 

The contribution from the wings of the PSF of the nu- 
cleus can be estimated by scali ng t he spectrum obtained 
from the nuclear readout trail (§ |4.3| ) by the fraction of the 
PSF falling within the chosen annulus. However, this scal- 
ing preserves the shape of the nuclear spectrum and does 
not account for the spectral flattening that occurs because 
of the energy dependence of the PSF. The deep calibration 
observation of the point source LMC X-l (obsid 1422) was 
used to model this effect. The readout trail image of LMC 
X-l was extracted along with spectra from various annuli 
surrounding the source. The channel-by-channel ratio of 
these spectra shows a linear rise with energy reflecting the 
PSF energy dependence. Applying this result to the nu- 
cleus of M81 modifies the power law form of the readout- 
trail spectrum to a function of the form E~ T (aE + b) with 
the constants a and b dependent on the annular region un- 
der study and the power law index T determined by the 
shape of the readout trail spectrum ( |4.3| ). Subsequently, 
all parameters for the nuclear contribution to the unre- 
solved emission are held constant while fitting the latter 
spectra. 

The contribution from unresolved sources is assumed 
to have the same shape as that of the resolved sources 
but with the model normalization left as a free parame- 
ter. The spectra of the resolved bulge sources, with the 
exception of the nucleus and the bright supersoft source, 
were added to obtain the total bulge source spectrum 



3 from archival Ultraviolet Imaging Telescope data available from http://archive.stsci.edu/uit/index.htm' 
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and then fitted to an absorbed power law. The result- 
ing column density is N 2 q = 7.7 and the photon index is 
T = 1.6 (x 2 = 156.6 for 176 dof). This slope is typical of 
individual intermediate-brightness sou rces for which reli- 
able spectral fits could be made (§ |3.3[ ). The total flux is 
9.2 x 10~ 13 ergs cm~ 2 s _1 corresponding to a luminosity 
of 1.43 x 10 39 ergs s _1 (1.64 x 10 39 ergs s -1 unabsorbed). 

The third component is assumed to be thermal because, 
when the data is fit with only the first two components, 
the bulk of the residual lies at low energies suggestive of a 
thermal contribution. Models with additional power law 
or thermal components were also tested. These compo- 
nents might represent a contribution from weak unresolved 
sources dissimilar in spectral shape to the resolved sources. 
None of these added components improved the fit statistic 
significantly. 

A typical spectrum and best-fit model components are 
displayed in Figure [ll]. For this annulus, the hard X-ray 
flux above ~1.5 keV is dominated by the nuclear contri- 
bution but the thermal and power law contributions are 
clearly present. 
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Fig. 11. — Spectrum of the unresolved emission from an annulus 
spanning 10 to 30" radius from the center of M81. Shown are the 
contributions from the nucleus, a power law representing unresolved 
point sources, and a thermal component. 

Figure |l2] shows that the temperature profile of the 
thermal component decreases with distance from the nu- 
cleus {(kT) = 0.38 ± 0.07 keV within 60" radius and 
(kT) = 0.26 ± 0.05 keV between 60" and 120", \ 2 = 7.2 
for 5 dof) and that the absorbing column density increases 
away from the nucleus (from N 2 o = 9.1 ± 2.6 to 21.2 ± 5.3 
on the same ranges, respectively, \ 2 — 13-2 for 5 dof). 
Assuming a homogeneous distribution of hot gas is re- 
sponsible for the thermal emission, its density and mass 
can be estimated based on the shape of the X-ray emit- 
ting region discussed above, namely a spherical bulge in 
the inner ~1 kpc surrounded by a region dominated by 
disk emission. The spectral fit parameters imply the num- 
ber density of the gas is n e ~ 0.01 cm -3 but rises to 
n e ~ 0.07 cm -3 in the innermost annulus and that the 
total mass of hot gas within the volume extending ~2' 
from the nucleus is ~ 7 x 10 6 M . This is a crude esti- 
mate because the thickness of the disk is unknown and was 



taken to be 1 kpc or roughly the thickness of the Galac- 
tic disk. The derived mass is a small fraction of the total 
mass of the bulge, - 3 x 10 10 M (Tenjes et al. 1998). 
The total thermal energy in the hot gas is ~7 x 10 54 erg 
or about 0.3% of the kinetic energy in bulge stellar motion 
(see Tenjes et al. 1998). 
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Fig. 12. — Radial variation of the best-fit temperature (top) and 
column density (bottom) for the thermal contribution to the unre- 
solved emission from the bulge. 



Page et al. (2002a, b) argued that the unresolved emis- 
sion is from collisionally excited, optically thin plasmas 
based the RGS data obtained by XMM-Newton. A model 
consisting of 3 mekal components (with kT = 0.18 ± 
0.04, 0.64±0.04, and 1.7±0.2 keV, respectively) and an ab- 
sorbed power-law (with T « 1.95, representing the nuclear 
contribution) provides a good fit to the RGS data. From 
the ratio of the forbidden line to the intercombination lines 
in the O VII triplets they rule out the alternative photo- 
ionization model at the 95% confidence level. The ratio of 
the resonance line to the other lines in the O VII triplet 
also put an upper limit of 10 9 cm -3 to the electron-number 
density of the line emitting gas. This result together with 
the luminosity of diffuse gas inferred from the Chandra 
data implies that the gas is not uniformly distributed and 
it fills only a small fraction of the available volume. 

Figure [13| displays the individual contributions to the 
X-ray luminosity of the unresolved excess emission as a 
function of the distance from the center of the galaxy. 
Both the thermal and power law contributions are rel- 
atively flat compared to the nuclear contribution whose 
spatial-dependence is dictated by the shape of the PSF. 
Also shown in Figure [l^ is the possible contribution from 
thermal emission at the center of the galaxy ded uced from 
analysis of the nuclear readout trail spectrum of §4.3. This 



emission is consistent with an extrapolation of the thermal 
and/or power law emission from the bulge and therefore 
is not necessarily intrinsic to the active nucleus. 
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Fig. 13. — Radial dependence of the X-ray luminosity of the un- 
resolved emission from the bulge. Shown are the contributions from 
thermal (dashed), power law (dot-dashed), and nuclear PSF (dotted) 
model components. The solid line represents the total luminosity. 
Also shown is the estimated thermal emission component at the 
nucleus obtained by adding a thermal model to the_model of the 
spectrum extracted from the nuclear readout trail (§4.3) 



The total absorbed luminosity in the power law com- 
ponent is 5.8 x 10 38 ergs s" 1 and 4.2 x 10 38 ergs s _1 
in the thermal component. The total is less than the 
1.6 x 10 39 ergs s _1 emitted by the discrete bulge sources. 

The expected contribution of unresolved discrete sources 
to the unresolved emission can be estimated by extrapolat- 
ing the observed luminosity function of the resolved bulge 
sources to lower luminosity. The bulge luminosity function 
is displayed in Figure |lj. Also shown are the luminosity 
functions of the subsets of bulge sources within 1' and be- 
yond 1' of the nucleus. The only systematic differences 
between the three functions is a flattening of the distri- 
bution at low luminosities for sources near the nucleus. 
This is caused by the loss of sensitivity in this region due 
to the large contribution of the nucleus to the underlying 
background. Tennant et al. (2001) pointed out the break 
in the background-subtracted bulge luminosity function at 
^4 x 10 37 ergs s _1 (^200 counts). The luminosity function 
shown here includes 13 more X-ray sources because of the 
larger ellipse used to define the bulge and is not corrected 
for background sources. Nevertheless, the break can still 
be seen. Another change in slope can be envisioned at 
about 60 counts. This change caused by the loss of sensi- 
tivity near the center of the bulge. The luminosity function 
is a power law, N(> C) = (194.1 ± 12)C-°- 50±0 - 02 between 
13 and 200 counts and N(> C) = (135.6 ± 5.0)C~ a37±a01 
between 13 and 60 counts. Extrapolating these two func- 
tions to lower luminosities implies unresolved sources of 
the type contributing to the resolved-source luminosity 
function account for only about 4 to 8% of the excess 
X-ray counts or an equivalent fraction of the luminosity 
(assuming the spectral shape is preserved). 



Fig. 14. — Luminosity function for bulge sources (solid). Also 
shown are the luminosity functions for sources lying within 1' of 
the nucleus (dashed) and outside (dotted) this radius. Note the 
lack of weak detected sources in the region near the nucleus due 
to the high contribution from the nucleus and from the X-ray ex- 
cess to the source background in this region. The break in the 
luminosity function at ~ 200 counts corresponds to a luminosity of 
~ 4 X 10 37 ergs s _1 . 

In summary, the spectrum of the unresolved emission 
probably has contributions from weak, unresolved sources 
and from diffuse hot gas. While spectral models show 
roughly equal contributions from both these components, 
extrapolation of the resolved-source luminosity function 
predicts only a small contribution from unresolved discrete 
sources. Perhaps another, distinct, population of weaker 
sources is present and contributes substantially to the un- 
resolved emission. These sources must be quite weak (and 
hence numerous) as no such population has been resolved 
in observations of the nearby galaxy M31 to a limiting lu- 
minosity of ~ 6 x 10 35 ergs s -1 (Shirey et al. 2001). The 
thermal component needed to fit the observed excess X- 
ray spectrum may also be from a collection of weak sources 
or from truly diffuse gas at kT ~ 0.4 keV. The mass and 
thermal energy content of any X-ray-emitting gas is small 
compared to the total mass and energy confined in the 
bulge. Unlike the distribution of far-UV light, both the 
unresolved emission and the resolved bulge point-source 
population trace the optical profile from old bulge stars 
and are confined by the same gravitational potential. In 
particular, there is no evidence of enhanced activity in the 
core such as heating by shocks as may explain the wispy 
Ha emission (Devereux et al. 1995, 1997). 

6. THE M81 DISK AND SPIRAL ARMS 

In contrast to the bulge, there is no measurable excess 
X-ray emission in the disk imaged by Chandra,. In partic- 
ular, there is no (unresolved) X-ray signature of the spiral 
arms which are seen clearly at other wavelengths and in 
X-rays in other galaxies such as M83 (Soria & Wu 2002). 

6.1. Spatial Distribution & Luminosity Function 

The spiral arms are perhaps most clearly seen in UV im- 
ages (compare, e.g., Reichen et al. 1994 or Hill et al. 1995) 
where the population of young stars trace the principal fea- 
tures of the spiral arms. The UV provides a high contrast 



17 



by not sampling the underlying disk component that is 
clearly seen in optical light. The archival near-UV image 
obtained by the Ultraviolet Imaging Telescope, with ap- 
proximately 3" resolution, was used to define the locations 
of the spiral arms. Distances between X-ray source posi- 
tions and the nearest spiral arm was measured. Sources 
within the 7/64 major-axis ellipse defining the bulge/disk 
interface (§|J) were excluded as were those sources exterior 
to the D25 isophote at 26/9. Source distances from spiral 
arms are displayed in Figure |l5| against the number of de- 
tected X-ray counts. Average distances in three luminosity 
ranges are also shown. Clearly, the average distance to a 
spiral arm decreases as the source luminosity increases. 
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Fig. 15. — Distribution of resolved X-ray sources in the disk. 
Shown is the distance from the nearest spiral arm against the 
(background-corrected) source counts. The line shows the average 
distance to spiral arms for sources in the 10-100, 100-1000, and 
>1000 count ranges (100 counts corresponds to an unabsorbed lu- 
minosity of ~ 2.4 X 10 37 ergs s" 1 ). 



72 X-ray sources located interior to the D25 isophote and 
exterior to the bulge/disk interface was binned into three 
groups with an equal number of sources in each group 
corresponding to sources within 0.30 kpc of a spiral arm, 
those 0.30 - 0.75 kpc, and those located >0.75 kpc from 
spiral arms. Note the strong dependence of the shape of 
the luminosity function on distance from spiral arms ex- 
hibited in Figure [l6|. Sources located near spiral arms are 
expected to be associated with the young population of 
stars recently formed following passage through the den- 
sity wave while those further from the spiral arms should 
be related to a relatively older population. Thus the ob- 
served X-ray luminosity functions for these populations 
suggests the younger population contains brighter objects 
and exhibits a luminosity function with a constant power 
law index (N(> C) = (75.7 ± 5.9)C~ a48±0 - 03 ). Those 
further from the spiral arms, in contrast, are dominated 
by weaker sources and their luminosity functions steepen 
above ^250 and ~100 counts, respectively (luminosities 
^6 x 10 37 ergs s _1 and ~2.4 x 10 37 ergs s _1 ). A power 
law fit to the full data range for sources between 0.30 
and 0.75 kpc gives N{> C) = (100.3 ± 10.8)C-°- 54±0 03 
and for those furthest from the spiral arms, N(> C) — 
(127.1 ±16.6)C-°- 69±0 04 . 

6.2. X-ray Spectra 

The spectra of the 44 resolved disk sources imaged on 
device S3, with the exception of SN 1993J and source X6, 
were added to obtain a representative disk source spec- 
trum to compare to the bulge source spectrum. An ab- 
sorbed power law with the addition of a thermal com- 
ponent was an improvement over an absorbed power law 
alone (A% 2 = 33 for 3 additional parameters). The re- 
sulting column density is A20 = 8.1, the photon index 
is r = 1.4, and the thermal component temperature and 
abundance are kT = 0.22 keV and Z = 0.2Z Q ( X 2 = 193 
for 172 dof). The thermal component accounts for ~3% 
of the total 0.3 - 8.0 keV flux. The power law slope is 
flatter than typical of individual sour ces f or which reliable 
spectral fits could be made (r = 1.5, §3.2) and of the com- 
posite spectrum of the resolved bulge sources (r = 1.6). 
The bulge spectrum also did not require a thermal com- 
ponent. The need for a thermal component suggests some 
of the weaker resolved disk sources may be X-ray-bright 
SNRs undetected in radio or in the survey of Matonick & 
Fesen (1997). Alternatively, a thermal contribution from 
black hole binaries in high soft states cannot be discounted. 
On the other hand, the thermal component contributes to 
the spectrum only at Fe L and at lower energies. A large 
range of column densities among the individual sources 
could also mimic this effect. 
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Fig. 16. — Luminosity function for disk sources. Shown are the 
luminosity functions for sources within 0.3 kpc of spiral arms (solid), 
sources between 0.3 and 0.75 kpc (dotted), and sources more than 
0.75 kpc away from spiral arms (dashed). 

The spatial dependence of the luminosity function for 
the disk sample is shown in Figure |l^. The sample of 



7. DISCUSSION 

7.1. The Discrete X-ray Source Population 

The X-ray spectra of the brightest sources in the M81 
field are predominantly moderately-absorbed power laws 
with photon index T ~ 1.5. A similar spectral shape repro- 
duces the combined bulge resolved-source spectrum while 
the combined disk source spectrum (for sources on S3) re- 
quires an additional weak thermal component. A power 
law is indicative of accreting X-ray binaries (XRBs) and a 
large population of bright power law sources is consistent 
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with surveys of our Galaxy and the Local Group where the 
bright X-ray source population is dominated by low-mass 
XRBs (e.g., Grimm et al. 2001). 

The scarcity of X-ray sources detected in the radio, or 
correlated with optically-selected SNRs, or exhibiting a 
strong thermal X-ray spectrum implies an insignificant 
number of the bright X-ray sources in M81 are SNRs. Su- 
pernova remnants are common in the Magellanic Clouds 
and in the solar neighborhood but are relatively less lu- 
minous and short-lived compared to typical XRBs. Per- 
haps some of the weaker disk sources are SNRs and ac- 
count for the thermal emission present in the combined 
disk source spectrum. Multi-wavelength observations of 
optically-identified extra-galactic SNRs have revealed that 
these sources are typically very weak X-ray and radio emit- 
ters (e.g. Pannuti et al. 2000, Lacey & Duric 2001, Pannuti 
et al. 2002). 

Four X-ray sources are coincident with known M81 glob- 
ular clusters 3 other X-ray sources are coincident with 
optically-bright objects with colors consistent with globu- 
lar clusters and two other X-ray sources have candidate op- 
tical counterparts that appear extended in Hubble images. 
Approximately 10% of Galactic globular clusters contain 
X-ray sources. In the case of M81, two surveys of the 
globular cluster population report a total of 139 globu- 
lar clusters including 98 within the Chandra field of view. 
Thus, 4% to 7% are coincident with X-ray sources. DiS- 
tefano et al. (2002) find 25% of confirmed clusters in their 
M31 field of view contain X-ray sources and 10% of all 
globular cluster candidates have X-ray sources. They also 
report that most of the luminous M31 X-ray sources are in 
globular clusters. In contrast, only one of the 11 brightest 
M81 X-ray sources, Einstein source X7, is coincident with 
a globular cluster. 

In general, there are remarkably few counterparts to 
the resolved X-ray sources identified in our assessment of 
the extensive literature and available archival images of 
M81. If the majority of the resolved sources are XRBs, 
then they have companion stars and accretion disks that 
may be detectable in optical light. The Hubble images 
of M81 approach a limiting magnitude of V ~ 27 mag or 
My ~ —0.5 mag. Thus only O and B main sequence or 
giant companions or highly-luminous accretion disks from 
XRBs in active states will be detectable by Hubble. Later- 
type companion stars will not appear in the optical data. 

Low-mass XRBs are common in galaxies because they 
are long-lived and slowly-evolving. The last encounter of 
M81 and its companion galaxy M82 occurred some 500 
Myr ago based on the study of the ages of young star clus- 
ters in M82 (de Grijs, O'Connell, & Gallagher 2001). If 
the onset of the last major star formation episode in the 
bulge of M81 was triggered by this encounter, then the 
most massive members of the current population of main- 
sequence stars have masses ^2.5 M© (Maeder & Meynet 
1988). If these constitute the population of companion 
stars in the currently-active XRBs, then the resolved bulge 
sources are mostly low-mass XRBs. They will not have de- 
tectable optical counterparts because they do not have O, 
B, or giant star companions (and because of the bright, 
amorphous optical background of the bulge). In contrast, 
on-going star formation along the spiral arms should pro- 
duce a population of high-mass XRBs with massive O and 



B star companions. This environment is, however, also the 
location of obscuring atomic and molecular gas. The few 
correlations with HII regions suggests some of the X-ray 
sources are located in star forming regions that may be 
populated by massive stars. The high percentage of Hub- 
ble potential counterparts identified in the disk relative to 
the bulge also suggests an abundance of early-type stars 
in the vicinity of the disk sources and, potentially, a pref- 
erence for high-mass XRB systems in this environment. 
This is consistent with the distribution of XRBs in our 
Galaxy where high-mass XRBs are concentrated towards 
the Galactic plane and along spiral arms while low-mass 
XRBs show a concentration towards the Galactic center 
(Grimm et al. 2001). 

7.2. The Brightest M81 Sources 

In-depth analysis of three of the 4 brightest sources in 
the M81 field was presented in §(|and of the third-brightest 
source in Swartz et al. (2002). Interestingly, all three of 
the brightest non-nuclear source are far from typical XRBs 
as seen in our Galaxy. SN 1993J is a supernova, Einstein 
source X6 is a rare ultra-luminous X-ray source with possi- 
ble optical and radio counterparts, and source number 132 
is an exceptionally-bright and hot supersoft source candi- 
date (Swartz et al. 2002). 

7.2.1. SN 1993J 

SN 1993J appears to be evolving as expected based on 
the standard CSM interaction model of Chevalier (1982; 
Fransson et al. 1996) though a complete picture incorpo- 
rating models of the exploding star and its pre-supernova 
environment awaits detailed numerical calculation. The 
X-ray properties of SN 1993J reported here provide an im- 
portant constraint on any future theoretical investigations 
because the X-ray light curve is declining steadily, even at 
~7 yr, whereas the most-detailed numerical simulations 
to date (Suzuki & Nomoto 1995) predicted the light curve 
would drop precipitously long ago unless the CSM were 
clumpy. A clumpy CSM would produce a varying light 
curve with episodes of high X-ray flux occuring whenever 
clumps are overtaken by the outgoing shock wave (Chugai 
1993). In this scenario, the CSM consists of a rarified wind 
embedded with relatively dense clouds and the X-rays em- 
anate from the shocked gas of the clouds with little or no 
reverse shock emission. This is not what is observed spec- 
troscopically. The spectrum of SN 1993J is best modeled 
with a combination of thermal emission from a reverse 
shock and a hard component from a forward shock. 

7.2.2. Einstein Source X6 

The multi-wavelength properties of Einstein source X6 
are intriguing. The X-ray spectrum of the source is best- 
fit with a disk blackbody model. In this model, the X- 
rays come from the inner portions of an accretion disk 
surrounding a compact object. The inferred mass of the 
central object is ~18 M assuming the innermost disk ra- 
dius derived from the model corresponds to the last stable 
Keplerian orbit of a non-spinning black hole. The X-ray- 
model-derived bolometric luminosity is near the Edding- 
ton limit for an object of this mass. The X-ray flux from 
X6 has been persistent throughout the >20 yr of observa- 
tion. 
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X6 is located within a 5"-diameter (~90 pc) SNR can- 
didate according to Matonick & Fesen (1997) based on a 
high [SII]/Ha ratio indicative of collisional excitation in 
the cooling region behind a SNR shock. There are no 
emission lines present in the X-ray spectrum of X6 and an 
optically-thin thermal plasma model is a notably poorer 
fit to the X-ray data. The X-ray morphology of X6 is that 
of a point source with no evidence for extension. Thus, no 
X-ray evidence, besides a steady flux, supports the conjec- 
ture that X6 is a SNR. 

A weak radio source is present at the location of X6. 
Synchrotron emission is observed at radio wavelengths 
from relativistic electrons accelerated in SNR shocks and 
in jets emanating from some (Galactic) XRBs. A radio (or 
optical) light curve of the source at the location of X6 has 
not yet been constructed. The radio source was present at 
3.6 cm in late 1994 but not seen in a 6 cm image taken in 
late 1999. Analysis of other radio images is in progress. If 
the radio source proves to be variable, then it is not from 
a SNR. The observed radio flux density of ~95 /iJy is typ- 
ical of, for example, Magellanic Cloud SNRs (Filipovic et 
al. 1998) after accounting for the disparate distances. In 
comparison, radio jets associated with Galactic XRBs are 
weaker except during extreme outbursts. 

There is also an optical point source coincident with X6. 
If associated with X6, the optical emission may either be 
from a moderately-massive, 09 - Bl, companion, which 
may be a Be star, or from the accretion disk itself but 
does not come from an extended source at Hubble resolu- 
tion. 

X6 can be compared to well-studied nearby XRBs. An 
example of a high-mass system with a massive compact 
accretor is Cyg X-l (e.g., van Paradijs 1995). The op- 
tical counterpart to Cyg X-l is a 09.7 supergiant with 
colors similar to those of the X6 counterpart but with the 
higher optical luminosity of a supergiant compared to a 
main-sequence star. The putative black hole in Cyg X-l 
exceeds 7 M and is most probably —16 M Q . Cyg X-l is 
a persistent X-ray source as is X6. It displays the charac- 
tersitic high-soft and low-hard states typical of black hole 
XRBs (Tanaka & Lewin 1995) and does not exceed an 
X-ray luminosity of ~ 2 x 10 38 ergs s _1 . Cyg X-l is radio- 
bright during its low-hard state with a flux of ~15 mJy or 
0.007 ytiJy if it were placed at the distance of M81. Scaling 
this value upward by the ratio of the X-ray luminosities 
of X6 to Cyg X-l in its low- hard state (—700) results in 
a radio flux density of only 6 /iJy which would not be 
detectable. 

An example of a low-mass system with a massive com- 
pact object and strong radio emission is the microquasar 
GRS 1915+105 (e.g., Mirabel & Rodriguez 1999). This 
system is a rapidly variable X-ray and radio transient 
reaching a peak X-ray luminosity of ~ 1.5 x 10 39 ergs s^ 1 in 
its high-soft state, comparable to X6, and an average lumi- 
nosity of ~ 3.7 x 10 38 ergs s _1 . High extinction along the 
line of sight to GRS 1915+105 obscures the optical coun- 
terpart and accretion disk. Near-infrared spectroscopy 
(Greiner et al. 2001), however, shows the companion to 
be a K - M main-sequence star and, along with the or- 
bital period, constrained the compact object mass to be 
14 + 4 M Q . GRS 1915+105 is a strong radio emitter. Scal- 
ing to the distance of M81 and to the X-ray luminosity of 



X6 (a factor of -13 when GRS 1915+105 is in its hard 
state) results in a radio flux density of about one-half the 
X6 value. 

Thus, while monitoring at many wavelengths is required 
before any definitive statment can be made, it is intriguing 
to consider X6 may be an X-ray- and radio-bright member 
of the class of microquasars (see Mirabel & Rodriguez 1999 
for a review) consisting of an accreting black hole with a 
radio-bright jet but with unusually-steady X-ray flux. 

7.2.3. The M81 Nucleus 

The X-ray properties of the nucleus of M81 arc diffi- 
cult to deduce from the present dataset because of severe 
pileup. A relatively weak spectrum extracted from the 
readout trail was analyzed and found to be a power law 
of photon index T = 1.98 ± 0.08, consistent with numer- 
ous previous X-ray studies. The presence of Fe Ka emis- 
sion could not be confirmed because of the lack of counts 
above —5 keV. Variability of the source also could not be 
assessed. However, it was shown, with the aid of the high 
angular resolution of the Chandra image, that the con- 
tribution to the nuclear spectrum from thermal emission 
is small or non-existent. Any thermal X-ray component 
present in the region is consistent with an extrapolation 
of the unresolved bulge emission observed surrounding the 
nucleus and extending over an —4 kpc diameter region. 

7.3. The M81 Bulge 

In addition to 53 X-ray sources resolved in the Chandra 
image, the bulge of M81 emits ~ 10 39 ergs s _1 in unre- 
solved emission. This is —12% of the total non-nuclear 
emission from the entire galaxy and is distributed over an 
— 2'-radius region centered on the nucleus. Both the re- 
solved sources and unresolved emission trace the optical 
light from the old population of bulge stars. 

If the unresolved emission is also produced by stellar 
systems, then they are systems distinct from the resolved 
sources because extrapolation of the luminosity function of 
the resolved sources contributes <10% of the unresolved 
emission. The possible X-ray-luminous stellar systems be- 
low the detection limit are massive OB stars, Be XRBs, 
CVs, RS CVn stars, and, at a lower luminosity, late-type 
stars. However, individual late-type stars have X-ray lu- 
minosities only of order a few 10 27 to a few 10 28 ergs s _1 , 
requiring some 10 12 stars to produce the unresolved emis- 
sion. 

Massive OB stars with colliding winds can be strong X- 
ray emitters but are rare. None are found in the bulge 
of M81 (Devereux, Ford, & Jacoby 1997). Be XRBs are 
young high-mass systems and the Be companion star is 
optically bright. They are therefore also unlikely to be 
abundant in the galactic bulge. 

CVs are short-period (typically <1 day) binaries con- 
sisting of a white dwarf and a late-type low-mass com- 
panion (Warner 1995). They arc numerous and are 
long-lived. The magnetic CVs, with a magnetic white 
dwarf, are known to have X-ray luminosities as high as 
— 10 32 ergs s _1 . The space density of magnetic CVs in the 
solar neighborhood is — 10~ 6 pc~ 3 (Warner 1995) while 
the stellar density is about 0.7 M pc -3 (Allen 1973). This 
implies the density of magnetic CVs is — 10~ 5 M Q _1 . If 
M81 has a similar space density of magnetic CVs, then 
there will be roughly 10 5 CV systems in the bulge of M81. 
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If about 10% are active (a rough estimate based on the 
properties of the local systems), then only ~ 10 36 ergs s _1 
of the unresolved bulge X-ray emission can come from 
CVs. 

RS CVn systems, composed of chromospherically ac- 
tive G or K stars with late-type main sequence or sub- 
giant companions, also have high X-ray luminosities. Typ- 
ical X-ray luminosities of RS CVn systems range from 
~ 10 29 ergs s" 1 to ~ 3 x 10 31 ergs s _1 (Rosner, Golub, & 
Vaiana 1985). Thus some 10 7 to 10 10 RS CVn systems are 
required to produce the unresolved bulge emission. If all 
stars in the bulge are ~1 M Q and half are in binary sys- 
tems, then there are ^10 10 binary systems in the bulge. 
Assuming about 20% of these systems become RS CVns 
and that G-K stars spend only a few percent of their life- 
times in their giant stage, an uncomfortably large fraction 
must currently be in an RS CVn phase. 

Individually, therefore, none of the these stellar sys- 
tems can readily account for the observed unresolved X-ray 
emission from the M81 bulge. If, instead, some portion of 
the unresolved bulge X-ray emission is from hot diffuse 
gas, as suggested by its spectral distinction from the sim- 
ple power law shape of the resolved sources, then only a 
small fraction (~0.02%) of the total bulge mass is needed 
to account for the observed emission. The X-ray emission, 
however, does not appear filamentary like the Ha emission 
does (Devereux et al. 1995). A filamentary morphology 
would be expected if the emission is from ionization by 
shocks. A source for producing shocks is also not obvious. 
Devereux et al. (1997) suggest shocks originating from nu- 
clear activity can account for the wispy "nuclear spiral" of 
Ha emission confined to the central but the unresolved 
X-ray emission is rather smoothly extended over a region 
of 2' radius. While ionizing radiation from hot evolved 
post-AGB stars may produce the observed UV excess in 
the core of M81 (O'Connell et al. 1992, Devereux et al. 
1995), these stars do not produce adequate ionizing radi- 
tion in the Chandra energy band to account for the X-ray 
emission (Binette et al. 1994). 

7.4. The M81 Disk and Spiral Arms 

One of the most spectacular features of M81 is its grand 
design spiral arm structure. The spiral arms trace the loca- 
tion of recent star-forming activity induced by the passage 
of spiral density waves. Applications of classical density 
wave models to M81 (e.g., Visser 1980) predict that mate- 
rial travels faster than the spiral pattern, entering an arm 
on the inside "upstream" edge. Stars forming at the spiral 
shock front travel at the local circular velocity of galactic 
rotation so that the youngest stars would be immediately 
downstream of the shock or toward the outside edge of the 
arm. The most massive stars are the quickest to evolve. 
They end their lives in core-collapse SN explosions leaving 
behind a neutron star or, perhaps, a black hole remant. 
Supernova explosions produce X-ray emitting SNRs and 
compact stars in binaries may become XRBs. Thus, the 
spiral arms are not only the site of star formation but also 
a stellar graveyard and the birthplace of X-ray sources. 

The brightest X-ray sources in the disk of M81 correlate 
spatially with the spiral arms. Accepting that the major- 
ity of the resolved sources are XRBs and that the X-ray 
flux is generally proportional to the mass accretion rate, 
then the brightest XRBs are young high-mass XRBs with 



high mass-transfer rates. The onset of mass-transfer in 
these systems, and hence of the X-ray-bright phase, can 
begin immediately following the formation of the compact 
object because of the strong stellar wind from the mas- 
sive companion. This is in contrast to the low-mass XRBs 
in which mass transfer begins only after the companion 
star evolves to a (sub)giant stage or when the binary or- 
bit has decayed sufficiently so that Roche lobe overflow 
can begin. Thus, the young high-mass systems become 
X-ray emitters while still within the spiral arm region of 
their origin. For this reason the luminosity distribution of 
the young XRBs in the spiral arms are expected to differ 
from the distribution of the older XRBs elsewhere in the 
galaxy. In particular, it will not show the characteristic 
luminosity break induced by aging of the XRB population 
as predicted by Wu (2001) and Wu et al. (2002a,b). 

Core-collapse supernova only come from stars more mas- 
sive than ~8-10 M and are X-ray bright SNRs only for 
a relatively short time. They, too, should be found prefer- 
entially near their place of origin, the spiral arms. 

8. CONTRIBUTIONS TO THE M81 X-RAY LUMINOSITY 
FUNCTION 

The observed M81 X-ray luminosity functions reported 
by Tennant et al. (2001) prompted Wu (2001) and Wu et 
al. (2002a, b) to consider the physical underpinnings that 
give rise to a cutoff in the luminosity function for the bulge 
sources and to the absence of this feature in the disk pop- 
ulation. Wu (2001) showed that the shape of the luminos- 
ity function is governed to first order simply by the birth 
rates and (X-ray active) lifespans of the XRBs that dom- 
inate the luminosity function and hence is a measure of 
the star formation history of the local environment and of 
galaxy evolution in the broader context. Further investiga- 
tion (Wu et al. 2002a,b) reveals that several complicating 
issues must be considered before this hypothesis can be rig- 
orously applied. Some of these issues have been addressed 
in the present work. 

The first issue is the presence of a population of SNRs. 
The onset and duration of the X-ray active phase of XRBs 
depends mainly on the donor star mass and its consequent 
evolutionary path while the X-ray luminosity depends on 
the accretion rate. This is fundamentally different than 
the X-ray evolution of SNRs. Here we have shown, how- 
ever, that SNRs are not important contributors to the to- 
tal X-ray source population in M81 with the exception of 
SN 1993J, the fourth-brightest source in the M81 field at 
the time of observation. 

Another issue is the occurrence of XRBs in globular clus- 
ters. If capture processes govern the formation of XRBs in 
globular clusters, as seems likely to account for the excess 
of XRBs in these environments, then the characteristic 
lifetime of the XRB is not just correlated with the nuclear 
or orbital evolution timcscalcs of the system but is also 
a function of the encounter frequency. However, only a 
few percent of the X-ray sources in M81 appear to be in 
globular clusters. As with Galactic globular cluster XRBs 
(Verbunt & van den Heuvel 1995), the impact on the lumi- 
nosity function is further minimized by the fact that the 
globular cluster XRBs in M81 are not among the brightest 
X-ray sources. Again, the exception is Einstein source X7, 
the fifth-brightest source in the field. 

A third factor with potential impact to the basic hypoth- 
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esis of Wu (2001) is the presence of XRBs with a nuclear- 
burning white dwarf accretor, i.e., members of the class of 
supersoft sources (e.g., Kahabka & van den Heuvel 1997). 
While the lifespans of these objects depend on the com- 
panion mass and mass-transfer rates as with other XRBs, 
only a narrow range of mass-transfer rates result in steady 
nuclear burning. Wu et al. (2002a, b) argue, therefore, a 
narrow X-ray luminosity range for this source population 
and a sharp decline in the number of sources with luminosi- 
ties above the Eddington limit for a Chandrasekar-mass 
accretor. Swartz et al. (2002) found 9 supersoft source 
candidates in the M81 field. Six of these are relatively 
weak sources with luminosities in a narrow range around 
~ 10 37 ergs s -1 . The three brightest candidates, how- 
ever, radiate at or above the Eddington limit depending 
on the adopted spectral model with one source (the third- 
brightest in the entire field) approaching ~ 10 39 ergs s _1 . 
Thus the most important contribution to the luminosity 
function from supersoft sources is at the high luminosity 
end and is dominated by one bright source. The effect 
of the weaker supersoft sources is obscured by the large 
number of other sources contributing at low luminosities. 




10 100 1000 

Source Counts (S) 



Fig. 17. — Observed X-ray luminosity function including all 123 
non-nuclear X-ray sources detected within the D25 isophote of M81 
(heavy solid line). Also shown are the luminosity functions of the 
SNRs (dotted), of X-rays sources spatially-coincident with globular 
clusters dashed), and of the supersoft sources (dot-dashed). The 
luminosity function with these three populations and the Ultra- 
Luminous X-ray source X6 omitted is shown as a thin solid line. 
The symbol (*) marks the division between bright sources for which 
spectral analysis has been reported in this work and of weak sources 
for which no spectral fits were made. 



The luminosity function for all non-nuclear sources de- 
tected within the D25 isophote of M81 is shown in Fig- 
ure [I?]. Also shown are the luminosity functions of the su- 
persoft sources, of X-ray sources spatially-coincident with 
globular clusters, and with SNRs. As shown in Figure [l7[ 
these three populations all have relatively flat power-law 
luminosity functions and they affect only the bright end of 
the overall luminosity distribution. XRBs, the dominant 
population of X-ray sources in M81, however, have a steep 
luminosity function and hence determine the overall shape 
of the luminosity functions, especially at the faint ends. 
The break at the luminosity of ~ 4 x 10 37 erg s _1 that 
we have found (see also Tennant et al. 2001) is therefore a 
characteristic imprint of the XRBs. We have argued that 
the formation of such a break is due to the age of a popu- 
lation of XRBs which were born at a star-burst episode in 
the recent past (Wu 2001, Wu et al. 2002a,b). This break 
is distinguishable from another possible break, expected 
to occur at ~ 2 x 10 38 erg s— 1, the Eddington luminos- 
ity of a 1.5 M© accreting object. The latter is attributed 
(Sarazin, Irwin, & Bregman 2001) to the presence of a 
population of neutron stars which accrete at rates close to 
the Eddington limit. Whether or not this break is visible 
in a given population depends on the relative proportion 
of neutron-star XRBs and black- hole XRBs. It also re- 
quires that the host galaxies (e.g. giant ellipticals) have a 
sufficiently large X-ray source population that the break 
becomes statistically significant. Nevertheless, we do see 
hints of this break in the luminosity function of the X- 
ray sources in M81, when we remove the SNRs, globular 
clusters XRBs and the supersoft sources. 

In summary, most of the factors complicating the simple 
birth-death model are unimportant for M81. Nevertheless, 
careful examination of the brightest sources is warranted 
because they have the largest influence on the luminosity 
function and yet are certainly not typical of the dominant 
class of X-ray sources, the XRBs. 

We thank L. Townslcy for applying her CTI-corrector 
algorithm to the Chandra data and for providing response 
matrices. We thank J. Davis for discussions and his inde- 
pendent analysis of the spectrum of the nucleus. We are 
grateful to N. Bartel and M. Beitenholtz for sharing their 
radio data. T. G. P. acknowledges a travel grant from 
the NRAO to reduce the radio data and he is very grate- 
ful to M. Rupen for guidance during the reduction process. 
K. W. thanks M. Weisskopf for providing support for visits 
to MSFC. Support for this research was provided in part 
by NASA/ Chandra grants GO0-1058X and AR2-3008X to 
D. A. S. 



REFERENCES 



Adler, D. S., k, Westphal, D. J. 1996, AJ, 111, 735 

Allen, C. W. 1973 Astrophysical Quantities (London: Athlone) 

Arnaud, K. 1996, in Astronomical Data Analysis Software and 

Systems V, ASP Conference Series 101, ed. G. Jacoby & J. Barnes 

(San Francisco:ASP) 17 
Bartel, N., et al. 1982, ApJ, 262, 556 
Bartel, N., et al. 2000, Sci., 287, 112 
Bash, F. N., U Kaufman, M. 1986, ApJ, 310, 621 
Bietenholz, M. F., Bartel, N., & Rupen, M. P. 2001, ApJ, 557, 770 
Binette, L., Magris, C. C, Stasinska, C, & Bruzual, A. G. 1994, 

A&A, 292, 13 
Bower, G. A., et al. 2000, BAAS, 197, 92.03 
Bradt, H. V. D., & McClintock, J. E. 1983, ARA&A, 21, 13 



Chandar, R., Ford, H. C, k, Tsvetanov, Z. 2001, AJ, 122, 1330 
Chevalier, R. A. 1982, ApJ, 258, 790 
Chevalier, R. A., & Blondin, J. M. 1995, ApJ, 444, 312 
Chugai, N.N. 1993, ApJ, 414, L101 
Cottrell, G. A. 1977, MNRAS, 178, 577 
Dalton, W. W., & Sarazin, C. L. 1995, ApJ, 440, 280 
de Grijs, R., O'Connell, R. W., & Gallagher, J. S. Ill 2001, AJ, 121, 
768 

de Jong, R. S. 1996, A&AS, 118, 557 

de Vaucouleurs G., de Vaucouleurs A., Corwin Jr. H. G., Buta R. 

J., Paturel G., &; Fouque P. 1991, Third Reference Catalogue of 

Bright Galaxies (NewYork: Springer- Verlag) 
Davidge, T. J., & Courteau, S. 1999, AJ, 117, 2781 



22 



David, L. P., Jones, C, & Forman, W. 1992, ApJ, 388, 82 
Davis, J. E. 2001, ApJ, 562, 575 

Devereux, N. A., Ford, H., & Jacoby, G. 1997, ApJ, 481, L71 
Devereux, N. A., Jacoby, G., & Ciardullo, R. 1995, AJ, 110, 1115 
Di Stefano, R., Kong, A. K. H., Garcia, M. R., Barmby, P., Greiner, 

J., Murray, S. S., & Primini, F. A. 2002, ApJ, 570, 618 
Elvis, M., & van Speybroeck, L., 1982, ApJ, 257, L51 
Fabbiano, G. 1988, ApJ, 325, 544 
Fabbiano, G. 1989, ARA&A, 27, 87 
Ferrarese, L., et al. 2000, ApJ, 529, 745 
Filipovic, M. D., et al. 1998, A&AS, 127, 119 

Fransson, C, Lundqvist, P., & Chevalier, R. A. 1996, ApJ, 461, 993 
Freedman, W. L., et al. 1994, ApJ, 427, 628 

Freeman, P. E., Kashyap, V., Rosner, R., & Lamb, D. Q. 2002, ApJS, 
138, 185 

Ghosh, K. K., Swartz, D. A., Tennant, A. F 

380, 251 
Goad, J. W. 1976, ApJS, 32, 89 
Greenwalt, B., Walterbos, R. A. M., Thilker, D 

1998, ApJ, 506, 135 
Greiner, J., Cuby, J. G., McCaughrean, M. J. 

& Mennickent, R. E. 2001, A&A, 373, L37 . 

Grim m, H , J ,., Gilfanov, M., & Sunyaev, R. 2001, preprint (astro 



Wu, K. 2001, A&A, 



Hoopes, C. G. 



Castro-Tirado, A. J. 



i Jacoby, G. H. 1996, AJ, 
Prestwich, A. H., Raley, 



ph/C 109239) 
Helfand, D. J 1984, PASP, 96, 913 
Hill, J. K., et al. 1992, ApJ, 395, L37 
Hill, J. K., et al. 1995, ApJ, 438, 181 

Ho, L. C, Filippcnko, A. V., & Sargent, W. L. 1996, ApJ, 462, 183 

Hodge, P. W., & Kennicutt, R. C, Jr. 1983, ApJ, 267, 563 

Houck, J. C, & Fransson, C. 1996, ApJ, 456, 811 

Immlcr, S., Aschcnbach, B., & Wang, Q. D. 2001, ApJ, 561, L107 

Immler, S., & Wang, Q. D. 2001, ApJ, 554, 202 

Ishisaki, Y., et al. 1996, PASJ, 48, 237 

Ivanova, G. R., 1992, MNRAS, 257, 119 

Iyomoto, N., & Makishima, K. 2001, MNRAS, 321, 767 

Kahabka, P., & van den Heuvel, E. P. J. 1997, ARA&A, 35, 69 

Kaufman, M., Bash, F. N., Kennicutt, R. C, Jr., & Hodge, P. W. 

1987, ApJ, 319, 61 
Kaufman, M., Bash, F. N., Hine, B., Rots, A. H., Elmegreen, D. M., 

& Hodge, P W. 1989, ApJ, 345, 674 
Kaufman, M., Bash, F. N., Crane, P. C, & 

112, 1021 

Kilgard, R. E., Kaaret, P Kr a „ aa M T 
T., & Zezas, A. 2002, astro-ph/020319C| 

Kong, X., et al. 2000, A J, 119, 2735 

Lacey, C. K. & Duric, N. 2001, ApJ, 560, 719 

Leahy, D. A., Darbro, W., Eisner, R. F., Weisskopf, M. C, Kahn, S., 
Sutherland, P. G., & Grindlay, J. E. 1983, ApJ, 266, 160 

Long, K. S., & van Speybroeck, L. P. 1983, in Accretion-Driven X- 
Ray Sources, ed. W. H. G. Lewin, & E. P. J. van den Heuvel 
(Cambridge: Cambridge Univ. Press) 117 

Ma, et al. 1998 Astron. J., 116, 516 

Maccacaro, T., Gioia, I. M., Wolter, A., Zamorani, G., & Stocke, J. 

T. 1988, ApJ, 326, 680 
Maeder, A., & Meynet, G. 1988, A&AS, 76, 411 
Makishima, K., et al. 2000, ApJ, 535, 632 

Maoz, D., Koratkar, A., Shields, J. C, Ho, L. C, Filippenko, A. V., 

& Sternberg, A. 1998, AJ, 116, 55 
Marcaide, J. M., et al. 1997, ApJ, 486, L31 
Marcaide, J. M., et al. 1993 IAUC Circ. 5820 

Matheson, T., Filippenko, A. V., Ho, L. C, Barth, A. J., & Leonard, 

D. C. 2000, AJ, 120, 1499 
Matonick, D. M., & Fesen, R. A. 1997, ApJS, 112, 49 
Mirabel, I. F., & Rodriguez, L. F. 1999, ARA&A, 37, 409 
Mitsuda, K., et al. 1984, PASJ, 36, 741 
Mizuno, T. 2000, Ph.D. thesis, Univ. Tokyo 
Nahar, S. N. 1999, ApJS, 120, 131 

Nandra, K., George, I. M., Mushotzky, R. F., Turner, T. J., & 

Yaqoob, T. 1997, ApJ, 477, 602 
Nomoto, K., Iwamoto, K., & Suzuki, T. 1995, Phys. Rep., 256, 173 
O'Connell, R. W., et al. 1992, ApJ, 395, L45 

Ohashi, T., et al. 1992, in Windows on Galaxies, eds G. Fabbiano, 
J. Gallagher, & A. Renzini (Dordrecht:Kluwer), 243 

Page, M. J., Breeveld, A. A., Soria, S. S., & Wu, K., 2002a, in New 
Visions of the X-ray Universe in the XMM-Newton and Chandra 
Era, in press 

Page, M. J., et al. 2002b, A&A, submitted 

Pannuti, T. G., Duric, N., Lacey, C. K., Goss, W. M., Hoopes, C. 

G., Walterbos, R. A. M. & Magnor, M. A. 2000, ApJ, 544, 780 
Pannuti, T. G., Duric, N., Lacey, C. K., Ferguson, A. M. N., Magnor, 

M. A. & Mendelowitz, C. 2002, ApJ, 565, 966 
Park, S., et al. 2002, ApJ, 567, 314 



Patat, F., Chugai, N., & Mazzali, P. A. 1995, A&A, 299, 715 
Peimbert, M., & Torres-Peimbert, S. 1981, ApJ, 245, 845 
Pclligrini, S., Cappi, M., Bassani, L., Malaguti, G., Palumbo, G. G. 

C, & 

Persic, M. 2000, A&A, 353, 447 

Perelmuter, J. -M., Brodie, J. P., & Huchra, J. P. 1995, AJ, 110, 620 
Perelmuter, J. -M., & Racine, R. 1995, AJ, 109, 1055 
Petit, H., Sivan, J. -P., & Karachentsev, I. D. 1988, A&AS, 74, 475 
Petre, R., Mushotzky, R. F., Serlemitsos, P. J., Jahoda, K., Marshall, 

F. E. 1993, ApJ, 518, 644 

Podsiadlowski, Ph., Hsu, J. J. L., Joss, P. C, & Ross, R. R. 1993, 
Nature, 364, 509 

Quataert, E., di Matteo, T., Narayan, R., & Ho, L. C. 1999, ApJ, 
525, L89 

Radecke, H. -D. 1997, A&A, 319, 18 

Reichen, M., Kaufman, M., Blecha, A., Golay, M., & Huguenin, D. 

1994, A&AS, 106, 523 
Rice, W. 1993, AJ, 105, 67 

Roberts, M. S. 1972, in IAU Symp. 44, External Galaxies and Quasi- 
Stellar Objects (Holland:Reidel), 12 
Roberts, W. W., Jr., & Hausman, M. A. 1984, ApJ, 277, 744 
Rosner, R., Golub, L., & Vaiana, G. S. 1985, ARA&A, 23, 413 
Rots, A. H., & Shane, W. W. 1975, A&A, 45, 25 
Sarazin, C. L., Irwin, J. A., & Bregman, J. N. 2001, ApJ, 556, 533 
Shara, M. M., Sandage, A., & Zurek, D. R. 1999, PASP, 111, 1367 
Shirey, R. et al. 2001, A&A, 365, L195 

Sholukhova, O. N., Fabrika, S. N., Vlasyuk, V. V., & Dodonov, S. 

N. 1998, AZh, 24, 507 
Soria, R., & Wu, K. 2002, A&A, 384, 99 

Stark, A. A., Gammie, C. F., Wilson, R. W., Bally, J., Linke, R. A., 

Heiles, C, & Hurwitz, M. 1992, ApJS, 79, 77 
Stetson, P. B. 1987, PASP, 99, 191 
Suzuki, T., & Nomoto, K. 1995, ApJ, 455, 658 

Swartz, D. A., Ghosh, K. K., Suleimanov, V., Tennant, A. F., Wu, 

K., 2002, ApJ, in press 
Tanaka, Y. & Lewin, W. H. G. 1995, in X-ray Binaries, ed. W. H. 

G. Lewin, J. van Paradijs, & E. P. J. van den Heuvel (Cambridge 
University Press, Cambridge) 126 

Tenjes, P., Haud, U., & Einasto J. 1998, A&A, 335, 449 

Tennant, A.F., Wu, K., Ghosh, K. K., Kolodziejczak, J. J., & Swartz, 

D. A. 2001, ApJ, 549, L43 

Terashima, Y., Iyomoto, N., Ho, L. C, Ptak, A. F. 2002, ApJS, 139,1 
Thielemann, F. -K., Nomoto, K., & Hashimoto, M. 1996, ApJ, 460, 
408 

Townsley, L. K, Broos, P. S., Garmire, G. P., & Nousek, J. A. 2000, 

ApJ, 534, L139 
Turner, T. J., & Pounds, K. A. 1989, MNRAS, 240, 833 
Uno, S., et al. 2002, ApJ, 565, 419 

Van Dyk, S. D., Weiler, K. W., Sramek, R. A., Rupen, M. P., & 

Panagia, N. 1994, ApJ, 432, L115 
van Paradijs, J. 1995, in X-ray Binaries, ed. W. H. G. Lewin, J. van 

Paradijs, & E. P. J. van den Heuvel (Cambridge University Press, 

Cambridge) 536 

van Paradijs, J. & McClintock, J. E. 1995, in X-ray Binaries, ed. W. 

H. G. Lewin, J. van Paradijs, & E. P. J. van den Heuvel (Cambridge 
University Press, Cambridge) 58 

Verbunt, F. & van den Heuvel, E. P. J. 1995, in X-ray Binaries, 
ed. W. H. G. Lewin, J. van Paradijs, & E. P. J. van den Heuvel 
(Cambridge University Press, Cambridge) 457 

Visser, H. C. D. 1980, A&A, 88, 159 

Wang, Q. D. 1999, ApJ, 517, 27 

Warner, B. 1995, Cataclysmic Variable Stars (New York: Cambridge 

University Press) 
Wheeler, J. C, & Filippenko, A. V. 1996, in Supernovae and 

Supernova Remnants, Proc. IAU Coll. 145, ed. R. McCray & Z. 

Wang (Cambridge: Cambridge Univ. Press), 241 
Woosley, S. E., Eastman, R. G., Weaver, T. A., & Pinto, P. A. 1994, 

ApJ, 429, 300 

Woosley, S. E., Pinto, P. A., & Hartmann, D. 1989, ApJ, 346, 395 
Wu, K. 2001, PASA, 18, 443 

Wu, K., Tennant, A. F., Swartz, D. A., Ghosh, K. K., & Hunstead, 
R. W. 2002a, in New Visions of the X-ray Universe in the XMM- 
Newton and Chandra Era, in press 

Wu, K., Tennant, A. F., Swartz, D. A., Ghosh, K. K., & Hunstead, 
R. W. 2002b, ApJ, submitted 

Xu, Y., Sutherland, P., McCray, R., & Ross, R. R. 1988, ApJ, 327, 
197 

Zickgraf, F. -J., & Humphreys, R. M. 1991, AJ, 102, 113 
Zickgraf, F. -J., Szeifert, Th., & Humphreys, R. M. 1996, A&A, 312, 
419 

Zimmermann, H. -U., et al. 1994, Nature, 367, 621 



Table 2 
M81 Discrete X-Ray Sources 
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9 


53 53.41 


RQ 

oy 


O^ P.Q A 

uo oy.4 


4.48 ± 0.94 


3.36 


— 


i3 


D 25 


0.82 ±0.17 




26 


9 


53 53.67 


RQ 

oy 


UO lo.U 


3.24 ±0.80 


2.87 


— 


i3 


D 25 


0.59 ±0.15 




27 


9 


53 57.47 


RQ 

oy 


UO OO.O 


16.92 ± 1.83 


7.49 


— 


i3 


D 25 


3.10 ±0.34 


H6, P15 


28 


9 


54 02.15 


RQ 

oy 


01 9R 7 
Ul ZO. / 


3.47 ±0.83 


3.02 


— 


i2 


D 25 


0.64 ±0.15 




29 


9 


54 06.73 


RQ 

oy 


08 A 1 7 
Uo 41. / 


11.37 ± 1.50 


5.95 


— 


i3 


d 


2.08 ± 0.28 




30 


9 


54 14.00 


RQ 

oy 


OP. 17 Q 

uo o ( .y 


12.96 ± 1.61 


6.71 


— 


i3 


d 


2.37 ±0.29 




31 


9 


54 14.92 


RQ 

oy 


OR 19 1 
UO 1Z.1 


2.75 ± 0.74 


2.84 


— 


i3 


d 


0.50 ±0.14 




32 


9 


54 21.30 


R52 
Do 


AA A 1 R 
44 41.0 


17.03 ± 1.84 


4.19 


— 


si 


D 25 


2.32 ± 0.25 




33 


9 


54 25.49 


R8 
Do 


/IR p.1 7 
40 Ol.i 


29.01 ±2.40 


6.75 


— 


si 


D 25 


3.94 ± 0.33 




34 


9 


54 26.42 


RS 
Do 


A'X A'X Pi 
4o 4O.0 


18.67 ± 1.93 


3.14 


— 


si 


D 25 


2.54 ±0.26 




35 


9 


54 27.84 


R8 
oo 


OO 11. o 


3.44 ± 0.83 


3.07 


— 


s2 


D 25 


0.63 ±0.15 




36 


9 


54 32.66 


oo 


A7 A A 1 

4 / 44. 1 


26.95 ±2.31 


5.77 


— 


si 


D 25 


3.66 ±0.31 




37 


9 


54 33.16 


R8 
oo 


p>9 9Q 
oz zy.u 


27.96 ±2.36 


9.51 


— 


s2 


D 25 


6.30 ± 1.10 


P18 


38 


9 


54 38.67 


oo 


p»9 49 Q 
oz 4z.y 


8.18 ± 1.28 


4.90 


— 


s2 


D 2h 


1.50 ±0.23 




39 


9 


54 39.23 


R8 
oo 


AP. zLQ A 

40 4y.4 


73.05 ±3.81 


12.55 


— 


si 


D 2h 


9.93 ±0.52 




40 


9 


54 41.82 


RS 
Oo 


p;r A7 (\ 

00 4 / .0 


11.91 ± 1.54 


6.51 


— 


s2 


D 2h 


2.18 ± 0.28 




41 


9 


54 41.99 


69 


02 43.7 


11.05 ± 1.48 


6.38 


>HST 


s3 


d 


1.50 ±0.20 


P21 


42 


9 


54 44.34 


68 


56 11.0 


5.61 ± 1.06 


4.44 


>19.81 


s2 


D 2 h 


1.03 ±0.19 




43 


9 


54 45.30 


68 


56 58.6 


58.15 ±3.40 


14.30 


>13.81 


s2 


D 25 


10.64 ± 0.62 


H8, P22, * 


44 


9 


54 46.79 


69 


05 12.6 


3.58 ±0.84 


3.07 


23.01 


s3 


d 


0.49 ±0.11 




45 


9 


54 47.18 


69 


01 01.4 


4.64 ± 0.96 


3.77 




s3 


d 


0.63 ±0.13 




46 


9 


54 51.49 


68 


51 43.5 


9.50 ±1.37 


5.31 




s2 


£>25 


1.74 ± 0.25 




47 


9 


54 53.96 


68 


54 55.0 


6.40 ± 1.13 


4.62 




s2 


D 25 


1.17 ±0.21 




48 


9 


54 55.15 


69 


04 20.3 


2.99 ±0.77 


3.40 




s3 


d 


0.41 ± 0.10 




49 


9 


54 55.60 


68 


51 59.6 


4.15 ±0.91 


2.85 




s2 


£>25 


0.76 ±0.17 




50 


9 


54 56.05 


69 


05 17.4 


3.43 ± 0.83 


3.18 


23.79 


s3 


d 


0.47 ±0.11 





24 



Table 2 — Continued 
M81 Discrete X-Ray Sources 







R.A. 

(J2000) 


Dec. 
(J2000) 


Count Rate 
(10- 4 s- 1 ) 


S/N 


mv 


CCD 


region 


L x 

(10 37 erg s" 1 ) 


Comment 


51 


9 




^7 ^Q 

o i .oy 


ftQ 

oy 


n9 a i i 

UZ "4:1.1 


32.91 ±2.56 


10.82 




s3 


d 


5.00 ± 0.60 




52 


9 


DO 


nn 1 1 

UU. 11 


ftQ 

oy 


n7 4^ 9 
Ui 40. z 


189.55 ±6.14 


25.64 


22.7(1) 


s3 


d 


27.00 ± 0.60 


X10, H10, P25 


53 


9 


DO 


nn 98 

UU.Zo 


ftQ 

oy 


n/i ^ft q 
U4 oo.y 


2.61 ±0.72 


2.86 





s3 


d 


0.36 ±0.10 




54 


9 


DD 


UU.oO 


ftQ 

oy 


ni 48 Q 
ui 4o.y 


2.80 ± 0.75 


2.90 





s3 


d 


0.38 ±0.10 




55 


9 


OO 


nn Ak 

UU.40 


ft8 

Oo 


^ft ^9 8 

00 oZ.o 


2.42 ±0.69 


2.91 





s2 


d 


0.44 ± 0.13 




56 


9 


00 


ni nn 

Ul.UU 


Oo 


99 i 

OO ZZ.l 


10.48 ± 1.44 


6.11 





s2 


£>25 


1.92 ± 0.26 


* 


57 


9 


00 


ni n^ 


ftQ 

oy 


n7 97 1 

U / Z ( . 1 


55.15 ±3.31 


13.88 





s3 


d 


22.00 ± 2.80 




58 


9 


00 


ni /in 
ui.4u 


(ik 

Oo 


9Q 7 

Oo zy. / 


13.66 ±1.65 


6.50 





s2 


£>25 


2.50 ± 0.30 




59 


9 


00 


ni 

U1.00 


ftQ 

oy 


i n AO Q 
1U 4z.o 


9.00 ± 1.34 


5.03 





s4 


d 


1.65 ± 0.24 


P26 


60 


9 


00 


n9 ^7 

UZ.O / 


dk 

Oo 


^ft 91 9 
OO Zl.Z 


12.73 ±1.59 


6.78 





s2 


£>25 


2.33 ± 0.29 


Hll, P27, * 


61 


9 


00 


n^ 4^ 

U0.4O 


fik 

Oo 


/1/1 99 8 
44 ZZ.O 


85.31 ±4.12 


12.94 





si 


D 25 


11.60 ±0.56 




62 


9 


00 


n^ ^9 

UO.OZ 


(ik 

Do 


^8 ^9 1 
Oo OZ.l 


3.60 ±0.85 


3.11 





s2 


d 


0.66 ±0.15 


P28, HII 


63 


9 


00 


n^ 

U0.O4 


(10 

oy 


n/i n^ 7 

U4 uo. / 


10.52 ± 1.45 


5.79 





s3 


d 


1.43 ± 0.20 




64 


9 


00 


n8 qi 
uo.y i 


(ik 

Oo 


^7 99 Q 

Oi zz.y 


2.89 ± 0.76 


3.03 


>19.69 


s2 


d 


0.53 ±0.14 




65 


9 


00 


nn 98 
uy.zo 


(Ik 
Oo 


Q^ ft 
Oo oo.o 


10.34 ± 1.43 


5.94 





s2 


£>25 


1.89 ± 0.26 




66 


9 


00 


nn 

uy.oo 


(10 

oy 


Q7 /1Q /I 
Ui 4o.4 


21.26 ±2.06 


8.43 


23.19 


s3 


d 


2.89 ± 0.28 


HII 


67 


9 


00 


nn 77 
uy. / / 


ftQ 

oy 


OA C\7 8 
U4 Ul.o 


18.93 ± 1.94 


8.29 


>17.29 


s3 


b 


3.50 ± 0.70 


X2, H13, P29, GC? 


68 


9 


00 


nn sn 
uy.ou 


dO 

oy 


08 /I 
Uo o0.4 


5.57 ± 1.05 


2.97 


21.45 


s4 


d 


1.02 ±0.19 




69 


9 


00 


1 n 90 
lu.zy 


(10 

oy 


n^ n9 /i 

UO UZ.4 


201.59 ±6.33 


27.36 


22.79 


s3 


b 


26.70 ± 1.00 


X3, H14, P31 


70 


9 


00 


l n 71 

1U. ( 1 


(10 

oy 


n8 A1 7 

Uo 4o. / 


13.25 ±1.62 


6.23 


22.86 


s4 


d 


2.42 ± 0.30 


P30, SNR 


71 


9 


00 


1 1 ^i 


(ik 

Do 


^7 47 Q 

o / 4 ( .y 


2.92 ±0.76 


2.91 


— 


s2 


d 


0.53 ±0.14 




72 


9 


00 


10 A A 


(10 

oy 


ni 91 ^ 

Ul Zl.O 


2.53 ±0.71 


2.83 


>HST 


s3 


d 


0.34 ±0.10 




73 


9 


00 


1/119 
14. 1Z 


(10 

oy 


1 9 Id 1 
1Z oO.l 


12.74 ± 1.59 


5.56 





s4 


d 


2.33 ± 0.29 


P32 


74 


9 


00 


1 /I £1 


(10 

oy 


nft 4n /i 

UO 4U.4 


3.94 ±0.89 


3.37 


24.03 


s3 


b 


0.54 ± 0.12 




75 


9 


00 


1 ^ 99 
10. ZZ 


ftQ 

oy 


UO oo.U 


7.40 ± 1.21 


4.83 


24.36 


s3 


b 


1.01 ± 0.16 




76 


9 


00 


10.00 


(ik 

Oo 


^A 97 ^ 
04 Z ( .0 


4.70 ± 0.97 


4.18 





s2 


d 


0.86 ±0.18 




77 


9 


00 


1 ^ QQ 

io. yy 


(ik 

Oo 


oi oy.o 


9.70 ± 1.39 


5.38 





s2 


£>25 


1.77 ±0.25 




78 


9 


00 


1 Q 7(1 

iy. / o 


(10 

oy 


n7 7 

U * oo. i 


6.08 ± 1.10 


3.90 


22.05 


s3 


d 


0.83 ±0.15 


HII, SNR 


79 


9 


00 


1 Q Q^ 

iy.yo 


ftQ 

oy 


n^ ^9 n 
Uo oz.u 


2.43 ± 0.70 


2.99 


>HST 


s3 


b 


0.33 ± 0.09 


Radio 


80 


9 


00 


91 8^ 
Zl.oO 


ftQ 

oy 


n^ /1 /i o 

Uo 44. y 


3.52 ±0.84 


3.34 


>HST 


s3 


b 


0.48 ± 0.11 




81 


9 


00 


91 87 
Zl.CS/ 


ftQ 

oy 


n^ 99 ^ 

UO ZZ.O 


37.95 ± 2.75 


11.53 


>HST 


s3 


b 


4.80 ± 0.50 




82 


9 


00 


91 QQ 

zi.yy 


(10 

oy 


fift Q7 ft 
UO oi .D 


35.61 ±2.66 


11.08 


17.14 


s3 


b 


5.10 ± 1.10 


H16, P33, GC? 


83 


9 


00 


99 if; 
ZZ. 10 


(10 

oy 


n^ i n ^ 

UO 1U.0 


178.55 ± 5.96 


25.84 


>HST 


s3 


b 


22.00 ± 0.90 


H15, P34 


84 


9 


00 


9Q 71 
Zo. / 1 


00 


^8 48 Q 
Oo 4o.y 


3.61 ±0.85 


3.85 





s2 


d 


0.66 ±0.16 




85 


9 


oo 


9/i ^n 


ftQ 
o y 


04 Q 

U4 oy.o 


6.06 ± 1.10 


4.46 


>HST 


s3 


b 


0.82 ±0.15 




86 


9 


oo 


9/1 QfS 


ftQ 
u y 


HQ ^7 Q 
uy o / .y 


155.20 ±5.55 


23.19 





s4 


d 


35.00 ± 3.00 


X4, H17, P35 


87 


9 


oo 


9/1 41 

Z4.4_L 


ftQ 
oy 


1 4 ^n 7 

14 OU. / 


12.24 ± 1.56 


4.80 


>15.82 


s4 


d 


2.24 ± 0.29 




88 


9 


oo 


9/1 77 
Z4. / / 


ftQ 
oy 


ni n 4 

Ul 10.4 


596.32 ± 10.89 


46.98 


>20.22 


s3 


d 


48.00 ± 2.00 


SN 1993J, H18, P36 


89 


9 


oo 


0F> *31 

ZU.Ol 


ftQ 
oy 


04 Q7 Q 

U4 OI .O 


12.55 ± 1.58 


6.52 


>HST 


s3 


b 


1.71 ± 0.21 




90 


9 


00 


9ft ^7 
ZO.O / 


ftQ 

oy 


n/i nn /i 

U4 UU.4 


3.82 ±0.87 


2.83 


>HST 


s3 


b 


0.52 ± 0.12 




91 


9 


55 


26.93 


69 


05 42.4 


6.29 ± 1.12 


4.18 


— 


s3 


b 


0.86 ±0.15 




92 


9 


55 


27.01 


69 


04 15.3 


16.24 ± 1.80 


7.45 


>HST 


s3 


b 


2.21 ± 0.24 




93 


9 


55 


27.28 


69 


02 48.0 


28.07 ± 2.36 


9.77 


>27.76 


s3 


b 


3.70 ± 0.50 




94 


9 


55 


27.85 


68 


49 52.9 


10.25 ± 1.43 


4.05 




si 


£>25 


1.39 ±0.19 




95 


9 


55 


28.03 


69 


04 07.9 


38.02 ± 2.75 


11.48 


>HST 


s3 


b 


5.70 ± 0.60 




96 


9 


55 


28.44 


69 


02 44.5 


12.02 ± 1.55 


6.04 


>HST 


s3 


b 


1.64 ± 0.21 


sss 


97 


9 


55 


28.82 


69 


06 12.9 


11.87 ± 1.54 


6.17 




s3 


b 


1.61 ± 0.21 


H20 


98 


9 


55 


29.20 


69 


03 21.1 


4.85 ±0.98 


4.11 


>HST 


s3 


b 


0.66 ±0.13 




99 


9 


55 


29.28 


69 


15 57.2 


9.92 ± 1.40 


3.21 




s4 


£>25 


1.82 ± 0.26 


H19 


100 


9 


55 


30.21 


69 


03 18.4 


12.64 ± 1.59 


6.47 


>HST 


s3 


b 


1.72 ± 0.22 





Table 2 — Continued 
M81 Discrete X-Ray Sources 







R.A. 

(J2000) 


Dec. 
(J2000) 


Count Rate 

(nr 4 s- 1 ) 


S/N 


mv 


CCD 


region 


L x 

(10 37 erg s" 1 ) 


Comment 


101 


9 


55 30.25 


(10 

oy 




6.55 ± 1.14 


4.48 


>HST 


s3 


b 


0.89 ±0.16 




102 


9 


55 31.38 


(10 

oy 


OA 10 ^ 
U4 iy.o 


46.44 ± 3.04 


12.62 


>HST 


s3 


b 


6.80 ±0.60 




103 


9 


55 32.61 


(10 

oy 


o^ i °. o 

UO lo.U 


2.99 ±0.77 


2.85 


>HST 


s3 


b 


0.41 ±0.10 




104 


9 


55 32.66 


(10 

oy 


UZ o±.*± 


3.22 ±0.80 


2.95 


22.87 


s3 


b 


0.44 ±0.11 




105 


9 


55 32.99 


(10 

oy 


UU oo.o 


1614 ± 18 


73.39 


24.13 


s3 


d 


270 ± 10 


X6, H21, P37 


106 


9 


55 33.17 


oy 


o°. ^ i 


5964 ± 34 


140.0 


15.42 


s3 


b 


3400 ± 233 


X5, H22, P38 


107 


9 


55 33.92 


(10 

oy 


o°. a°. °. 


13.92 ± 1.66 


6.21 


>HST 


s3 


b 


1.89 ±0.23 




108 


9 


55 34.12 


(10 

oy 


07 1 °. 1 
Ui lo.l 


7.75 ± 1.24 


4.73 


— 


s3 


d 


1.05 ±0.17 




109 


9 


55 34.32 


(10 

oy 


uo ou.y 


45.90 ± 3.02 


9.31 


>HST 


s3 


b 


4.20 ±0.40 




110 


9 


55 34.56 


(10 

oy 


uo oy.u 


8.64 ± 1.31 


4.35 


19.52 


s3 


b 


1.18 ±0.18 




111 


9 


55 34.62 


(10 

oy 


09 n 
UZ ou.u 


16.79 ±1.83 


7.58 


>HST 


s3 


b 


2.28 ±0.25 




112 


9 


55 34.65 


(10 

oy 


01 a 

Uo 01.4 


52.38 ± 3.23 


8.47 


22.73 


s3 


b 


5.80 ±0.60 




113 


9 


55 34.71 


(10 

oy 


OA ^ Q 

U4 oo.y 


38.03 ± 2.75 


11.84 


>HST 


s3 


b 


4.30 ±0.50 




114 


9 


55 34.81 


(10 

oy 


Uo lo.O 


16.21 ± 1.79 


7.23 


20.78 


s3 


b 


2.20 ±0.24 




115 


9 


55 34.90 


(10 

oy 


OA 07 O 

U4 u < .y 


19.51 ± 1.97 


7.17 


>HST 


s3 


b 


3.10 ±0.50 




116 


9 


55 34.98 


dO 

oy 


no AO o 
Uo 4Z.O 


118.14 ±4.85 


20.38 


>HST 


s3 


b 


19.50 ± 1.00 




117 


9 


55 35.28 


(ik 

Do 


OO 1U.0 


31.02 ±2.48 


10.59 





s2 


d 


6.00 ± 1.00 


H24 


118 


9 


55 35.29 


(10 

oy 


o°. i ^ q 
Uo io. y 


56.97 ±3.37 


13.90 


>HST 


s3 


b 


7.40 ± 0.50 


H23 


119 


9 


55 35.40 


(10 

oy 


O^ p.7 7 
UO Oi.i 


4.55 ±0.95 


3.64 





s3 


b 


0.62 ±0.13 




120 


9 


55 35.56 


(10 

oy 


o°. °. 

Uo 04.O 


46.87 ± 3.05 


10.55 


>HST 


s3 


b 


3.60 ± 0.60 




121 


9 


55 35.72 


(10 

oy 


0(1 "37 7 
UO ol.l 


7.65 ±1.23 


3.86 


>20.66 


s3 


d 


1.04 ±0.17 




122 


9 


55 36.29 


(10 

oy 


09 A A k 
UZ 44. 


13.49 ± 1.64 


6.76 


>HST 


s3 


b 


1.83 ±0.22 




123 


9 


55 36.45 


(10 

oy 


09 AO ^ 
UZ 4U.0 


7.69 ± 1.24 


4.67 


>HST 


s3 


b 


1.05 ±0.17 




124 


9 


55 36.75 


(10 

oy 


0(1 °.°v 9 
UO oo.Z 


13.82 ± 1.66 


6.90 


>19.83 


s3 


d 


1.88 ±0.23 


HII, SNR 


125 


9 


55 36.87 


(Ik 
Do 


OO OO.Z 


3.03 ± 0.78 


3.40 





s2 


d 


0.55 ±0.14 




126 


9 


55 37.05 


(10 

oy 


U4 oo.o 


28.54 ± 2.38 


10.34 


>HST 


s3 


b 


6.40 ± 1.00 




127 


9 


55 37.28 


(10 

oy 


09 07 °. 
UZ U / .o 


3.50 ±0.83 


3.11 


18.41 


s3 


b 


0.48 ±0.11 


GC? 


128 


9 


55 37.60 


(10 

oy 


OA p.7 7 
U4 O/.i 


10.60 ±1.45 


6.08 


23.10 


s3 


b 


1.44 ±0.20 




129 


9 


55 37.66 


(10 

oy 


o°. i a 9 

Uo io.z 


9.06 ±1.34 


5.43 


>HST 


s3 


b 


1.23 ±0.18 


sss 


130 


9 


55 38.62 


dk 

Oo 


AO 99 Q 
4y zz.y 


10.98 ± 1.48 


4.09 





si 


£>25 


1.49 ±0.20 




131 


9 


55 40.69 


(10 

oy 


01 0^. O 
ui uo.u 


3.20 ±0.80 


3.21 


>HST 


s3 


b 


0.43 ±0.11 




132 


9 


55 42.21 


(10 

oy 


0°. Id °. 
Uo OO.O 


800.98 ± 12.62 


53.88 


21.67 


s3 


b 


30.00 ± 1.00 


SSS, H25 


133 


9 


55 42.86 


(10 

oy 


O*} 07 (1 
Uo U < .0 


21.84 ± 2.08 


8.82 


— 


s3 


b 


3.60 ± 1.70 




134 


9 


55 43.17 


(10 

oy 


OA A ^ O 
U4 40. U 


5.91 ± 1.08 


4.21 


>HST 


s3 


b 


0.80 ±0.15 




135 


9 


55 43.34 


o y 


OA 9^ 9 

U4: ZO.Z 


2.92 ±0.76 


3.33 


>HST 


s3 


b 


0.40 ±0.10 




136 


9 


55 43.76 


Rk 
oo 


PQ 04. & 

oy u4.o 


38.18 ±2.75 


11.44 





s3 


d 


4.00 ±0.40 


H27 


137 


9 


55 44.63 


(10 
oy 


1 O^i 9 
1U uo.z 


14.50 ±1.70 


6.80 


>20.77 


s4 


d 


2.65 ±0.31 




138 


9 


55 44.71 


(10 
oy 


0^1 ^4 ^ 

UO 0*1.0 


12.12 ± 1.55 


6.17 


22.77 


s3 


d 


1.65 ±0.21 




139 


9 


55 45.91 


o y 


0*} 00 A 

UO UU.4: 


14.70 ± 1.71 


7.13 


>17.60 


s3 


b 


2.00 ±0.23 




140 


9 


55 46.16 


dk 


co /in 7 

Oo 4U. / 


3.87 ±0.88 


3.27 


— 


s2 


d 


0.71 ±0.16 




141 


9 


55 47.05 


69 


05 51.1 


66.39 ± 3.63 


15.28 


18.63 


s3 


d 


11.80 ±1.30 


H28, HII, GC 


142 


9 


55 47.96 


68 


59 28.2 


6.11 ± 1.10 


4.52 





s3 


d 


0.83 ±0.15 


SSS, HII 


143 


9 


55 48.19 


68 


59 15.1 


4.22 ±0.92 


3.88 


>21.00 


s3 


d 


0.57 ±0.12 


SSS 


144 


9 


55 49.41 


68 


58 36.3 


52.30 ±3.22 


13.31 




s3 


d 


14.50 ±1.50 


H29, P40, SNR 


145 


9 


55 49.52 


69 


08 12.0 


78.12 ±3.94 


16.74 




s4 


d 


14.80 ± 1.60 


H30, P42 


146 


9 


55 49.87 


69 


05 32.0 


428.91 ± 9.23 


39.54 


20.77 


s3 


d 


59.80 ± 1.60 


X7, H31, P41, GC 


147 


9 


55 51.54 


69 


04 10.5 


2.61 ±0.72 


2.86 




s3 


b 


0.36 ±0.10 




148 


9 


55 51.58 


69 


07 43.1 


3.27 ±0.81 


3.26 


18.81 


s4 


d 


0.60 ±0.15 


HII, GC 


149 


9 


55 53.13 


69 


05 20.1 


35.05 ± 2.64 


11.05 


22.20 


s3 


d 


1.00 ±0.00 


SSS 


150 


9 


55 53.31 


69 


02 06.5 


14.22 ± 1.68 


7.02 




s3 


b 


1.93 ±0.23 





Table 2 — Continued 



M81 Discrete X-Ray Sources 






R.A. 




Dec. 


Count Rate 


S/N 


mv 


CCD 


region 


Lx 


Comment 






(J2000) 


(J2000) 


(10" 4 s" 1 ) 










(10 37 erg s _1 ) 




151 


9 


EE 
OO 


Do. Oo 


oy 


U4 


^/i s 


7.00 ± 1.18 


4.60 


24.28 


s3 


d 


0.95 ± 0.16 




152 


9 


EE 
00 


tx p: ni 
00. Ul 


oy 


09 

uz 


00. y 


4.30 ± 0.93 


3.77 




s3 


b 


0.59 ± 0.13 




153 


9 


EE 
00 


OO.O ( 


Do 


Oo 


Oo.O 


2.49 ± 0.70 


2.86 




s3 


d 


0.34 ±0.10 




154 


9 


r.r. 
00 


^ 7Q 
00. i y 


oy 


1 n 

1U 


OK ^ 
Uo.O 


9.00 ± 1.34 


5.31 




s4 


d 


1.65 ± 0.24 




155 


9 


EE 
00 


00. 1 1 


oy 


0'} 
Uo 


1 9 9 


6.46 ± 1.13 


4.87 




s3 


b 


0.88 ± 0.15 


SSS, Radio 


156 


9 


r.r. 
00 


^« 91 
00. Zl 


oy 


UO 


1 /I 7 
14. ( 


9.16 ± 1.35 


5.56 


23.56 


s3 


d 


1.25 ± 0.18 




157 


9 


KC 
00 


7/1 
OD. i 4 


oy 


Oft 
Uo 


H9 
UZ.O 


5.84 ± 1.08 


4.25 


20.5(1) 


s4 


d 


1.07 ±0.20 


Extended? 


158 


9 


r.r. 
00 


^8 fil 
Oo.Ol 


(10 

oy 


r » 

uo 


9fi 9 
zu.z 


64.07 ± 3.57 


15.11 


26.37 


s3 


d 


16.00 ± 1.70 


H32, GC 


159 


9 


r.r. 
00 


oy . 10 


oy 


or; 
uu 


1 7 /I 
1 f .4 


20.14 ± 2.00 


8.19 


25.28 


s3 


d 


2.10 ± 0.50 




160 


9 


01) 


ni Q7 
ui.y 1 


Do 


00 


oy.o 


44.41 ± 2.97 


12.81 




s3 


d 


2.00 ± 0.30 


H33, P43 


161 


9 


K« 
01) 


09 fiQ 

uz.oy 


Do 


oy 


^ 9 
00. z 


40.46 ± 2.84 


12.24 




s3 


d 


3.40 ± 1.10 


H34 


162 


9 


01) 


09 78 
UZ. i 


Do 


00 


1 1 

44. U 


12.85 ± 1.60 


6.88 




s3 


d 


1.75 ± 0.22 




163 


9 


01) 


0*3 1 ^ 
Uo. 10 


oy 


09 
UZ 


10. 


4.81 ± 0.98 


4.06 




s3 


d 


0.65 ± 0.13 




164 


9 


01) 


0*} 9Q 
uo.zy 


oy 


1 
Ul 


07 *3 
U/.o 


3.77 ± 0.87 


3.64 




s3 


d 


0.51 ± 0.12 




165 


9 


01) 


0/1 ^fi 
U4.o0 


oy 


1 1 

1 1 


P.Q 7 

oy. ( 


5.81 ± 1.07 


3.43 




s4 


d 


1.06 ± 0.20 




166 


9 


01) 


0/1 P.Q 

u4.oy 


Do 


00 


9 

oy.z 


2.90 ± 0.76 


3.24 




s3 


d 


0.39 ± 0.10 




167 


9 


56 


04.93 


69 


03 


43.7 


2.61 ± 0.72 


2.90 


>19.09 


s3 


d 


0.36 ± 0.10 


HII 


168 


9 


56 


06.07 


68 


59 


40.7 


5.66 ± 1.06 


4.54 




s3 


d 


0.77 ±0.14 


HII 


169 


9 


56 


06.09 


69 


08 


33.5 


6.72 ±1.16 


4.67 





s4 


d 


1.23 ±0.21 




170 


9 


56 


07.84 


69 


03 


25.2 


15.10 ±1.73 


7.08 




s3 


d 


2.05 ±0.24 




171 


9 


56 


09.05 


69 


01 


06.4 


110.56 ±4.69 


20.30 




s3 


d 


4.00 ±0.50 


SSS, H36, P44 


172 


9 


56 


09.48 


69 


12 


49.5 


22.46 ± 2.11 


8.32 




s4 


,025 


5.10 ±0.80 


H35, P45 


173 


9 


56 


13.74 


69 


06 


30.6 


21.10 ± 2.05 


8.07 


24.75 


s3 


d 


2.87 ±0.28 


H37, P46 


174 


9 


56 


14.21 


69 


02 


24.3 


4.02 ±0.89 


3.58 


M9.05 


s3 


d 


0.55 ±0.12 


SSS, HII 


175 


9 


56 


14.42 


69 


02 


47.9 


13.66 ± 1.65 


6.81 




s3 


d 


1.86 ±0.22 




176 


9 


56 


14.85 


69 


03 


37.7 


3.48 ± 0.83 


3.36 




s3 


d 


0.47 ±0.11 


HII 


177 


9 


56 


27.46 


69 


10 


12.1 


22.57 ±2.12 


8.46 




s4 


D 25 


3.10 ±0.50 





TABLE 3 
M81 Bright X-Ray Sources 





N H 
(10 2() cm- 2 ) 


Spectral 
Parameter" 


xVdof 


Lx 

(10 37 erg s _1 ) 


Pks 


37 


27.2 ± 19.0 


1.63 ±0.60 


2.03/5 


6.3 ±1.1 


74.50 


51 


14.9 ± 10.0 


1.76 ±0.42 


5.36/8 


5.0 ±0.6 


81.80 


52 


13.6 ±4.0 


1.59 ±0.14 


38.87/38 


27.0 ±0.6 


65.10 


57 


99.1 ±40.0 


1.90 ±0.51 


8.00/13 


22.0 ±2.8 


12.10 


67 


4.0 ±0.0 


1.31 ±0.28 


4.42/4 


3.5 ±0.7 


1.89 


69 


11.8 ±3.0 


1.47 ±0.11 


58.11/43 


26.7 ± 1.0 


58.50 


81 


13.5 ±7.0 


2.32 ±0.35 


14.01/7 


4.8 ±0.5 


45.90 


82 


10.0 ± 10.0 


1.43 ±0.55 


3.65/8 


5.1 ± 1.1 


53.40 


83 


14.4 ±4.0 


1.79 ±0.13 


60.26/35 


22.0 ±0.9 


92.60 


86 


4.0 ±0.0 


1.08 ±0.13 


42.43/35 


35.0 ±3.0 


4.03 


93 


15.6 ±7.0 


2.03 ± 0.45 


1.98/4 


3.7 ±0.5 


59.80 


95 


4.0 ±0.0 


1.61 ±0.17 


14.45/11 


5.7 ±0.6 


17.60 


102 


4.0 ±0.0 


1.80 ±0.24 


12.14/11 


6.8 ±0.6 


82.00 


109 


9.7 ±5.0 


2.20 ±0.33 


8.81/7 


4.2 ±0.4 


70.30 


112 


4.0 ±0.0 


1.17±0.16 


13.81/8 


5.8 ±0.6 


70.30 


113 


14.2 ± 10.0 


2.24 ±0.50 


6.54/8 


4.3 ±0.5 


9.76 


115 


7.7 ±7.0 


1.26 ±0.85 


2.22/4 


3.1 ±0.5 


25.40 


116 


8.7 ±3.0 


1.44 ±0.12 


32.71/30 


19.5 ±1.0 


68.70 


117 


16.0 ± 13.0 


1.92 ±0.43 


6.49/6 


6.0 ± 1.0 


58.20 


118 


5.2 ±3.0 


1.53 ±0.25 


13.28/12 


7.4 ±0.5 


99.10 


120 


11.5 ±5.0 


2.57 ±0.50 


8.00/5 


3.6 ±0.6 


98.40 


126 


4.0 ±0.0 


1.71 ±0.15 


30.69/14 


6.4 ± 1.0 


56.30 


133 


24.5 ± 20.0 


2.15 ±0.80 


1.16/6 


3.6 ± 1.7 


21.00 


136 


10.3 ±5.0 


1.95 ±0.47 


7.68/8 


4.0 ±0.4 


91.00 


141 


12.3 ± 7.0 


1.22 ±0.20 


8.84/15 


11.8 ± 1.3 


22.90 


144 


13.6 ±5.0 


2.61 ±0.39 


12.15/8 


14.5 ±1.5 


0.61 


145 


4.0 ±0.0 


1.41 ±0.18 


18.82/22 


14.8 ±1.6 


72.40 


146 


11.5 ±2.0 


1.37 ±0.08 


72.49/81 


59.8 ±1.6 


5.31 


158 


28.7 ± 11.0 


3.61 ±0.96 


19.33/13 


16.0 ±1.7 


35.50 


159 


17.4 ± 17.0 


2.13 ± 1.10 


6.88/3 


2.1 ±0.5 


38.10 


160 


4.0 ±0.0 


0.20 ±0.02 


11.89/7 


2.0 ±0.3 


37.70 


161 


38.3 ± 15.0 


0.57 ±0.14 


10.64/9 


3.4 ± 1.1 


71.70 


172 


10.0 ±6.3 


1.21 ±0.37 


4.37/7 


5.1 ±0.8 


2.58 



"Power law photon index, T, except for source numbers 160 and 161 
where spectral parameter is blackbody temperature, kT, in keV. 



